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Among  a number  of  modern  surface  sensitive  analytical  techniques. 
Auger  electron  spectroscopy  (AES)  has  been  extensively  applied  for  glass 
surface  research  because  of  its  superior  characteristics.  AES  has  been 
generally  considered  as  a "non-destructive"  technique;  however,  a number 
of  reports  have  shown  that  the  incident  electron  beam  can  alter  the  sur- 
face composition  of  the  solid  being  analyzed.  Compositional  changes  in 
glass  have  been  particularly  severe  in  glasses  containing  alkali  ions 
and  it  has  been  manifested  by  the  time-dependent  alkali  Auger  signal  loss 
during  electron  bombardment. 

The  present  work  describes  fundamental  mechanisms  of  alkali  Auger 
signal  loss  from  glass  surfaces  during  electron  bombardment.  It  has  been 
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shown  that  this  is  the  result  of  electric  field  enhanced  diffusion  and 
electron  stimulated  desorption.  Electric  field  enhanced  diffusion  was 
demonstrated  in  soda-silica  glass  thin  films  deposited  on  stainless  steel 
substrates.  The  amount  of  charge  trapping  was  controlled  by  the  experi- 
mental conditions  and  the  electric  field  strength  in  glass  was  evaluated 
to  be  approximately  10^  V/cm.  Electron  stimulated  desorption  was  studied 
from  soda-silica  glass  fibers  fractured  in  situ.  It  was  found  that  sodium 
and  oxygen  were  desorbed  during  electron  bombardment.  Electron  stimulated 

_3 

desorption  became  important  with  beam  current  densities  greater  than  10 
2 

A/ cm  , and  the  field  enhanced  diffusion  limited  further  increase  in  desorp- 

-2  2 

tion  of  sodium  with  a current  density  greater  than  5 x 10  A/cm  . The 

-19  2 

desorption  cross  section  for  sodium  was  found  to  be  about  10  cm  for 
the  primary  electron  energy  of  0.5  to  4.0  keV.  The  desorption  mechanism 
was  also  studied.  A Menzel-Gomer-Redhead  model  was  not  adequate  to  explain 
sodium  desorption  from  glass,  but  a core  hole  Auger  decay  model  was  the 
most  likely  mechanism.  Thereby  the  cross  section  was  correlated  with 
this  mechanism  and  reasonable  agreement  with  experimental  results  was  made. 

A comprehensive  theory  for  the  alkali  Auger  signal  loss  has  been 
developed.  This  was  applied  to  low  temperature  measurements  and  the  ori- 
gin of  the  incubation  time  (during  which  the  alkali  Auger  signal  was  nearly 
consistent)  became  apparent. 

AES  has  been  applied  to  the  investigation  of  microstructural  charac- 
terization of  lithia-silica  glass-ceramics.  The  rates  of  lithium  Auger 
signal  loss  were  influenced  by  the  presence  of  a crystalline  phase  and 
were  correlated  with  the  volume  percent  of  crystallization. 
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CHAPTER  I 
INTRODUCTION 

Overview 

The  variation  of  glasses  used  in  industrial  and  laboratory  appli- 
cations is  extensive.  For  many  uses  treatments  such  as  polishing, 
etching,  chemical  cleaning  and  heat  treatment  are  performed  to  control 
properties  related  to  reflectivity,  hardness,  chemical  durability  or 
adhesion.  For  these  uses,  knowledge  of  the  chemical  surface  composition 
would  help  to  make  more  reliable  predictions  of  properties.  Unfortunately, 
numerous  surface  phenomena  and  characteristics  are  not  yet  explained  due  to 
the  difficulties  of  these  measurements.  Classical  techniques  can  provide 
information  about  the  outer  few  thousand  Angstroms  of  glass;  however, 
these  techniques  are  generally  slow  and  have  poor  depth  resolution. 

In  recent  years  modern  surface  sensitive  analytical  techniques,  such 
as  Auger  Electron  Spectroscopy  (AES)  [1],  Electron  Spectroscopy  for 
Chemical  Analysis  (ESCA)  [2],  Secondary  Ion  Mass  Spectroscopy  (SIMS)  [3], 
Ion  Scattering  Spectroscopy  (ISS)  [4],  Sputter  Induced  Photon  Spectroscopy 
(SIPS)  [5]  and  Rutherford  Backscattering  Spectroscopy  (RBS)  [6]  have  been 
extensively  developed.  A common  feature  of  these  techniques  is  that  the 
sample  is  brought  into  a vacuum  chamber  where  it  is  bombarded  with  a beam 
of  energetic  particles  such  as  electrons,  ions  or  photons  and  the  inter- 
action of  these  beams  with  the  sample  produces  secondary  particles  which 
can  be  analyzed  to  obtain  information  from  a thin  surface  layer.  All  of 


1 


2 


these  techniques  except  RBS  are  sensitive  to  the  outermost  atomic  layers 
of  a solid ; therefore  any  change  in  the  concentration  or  chemical  state 
of  surface  species  will  be  immediately  evident. 

The  technique  of  Auger  Electron  Spectroscopy  is  used  extensively 
for  glass  surface  research  because  of  its  good  spatial  resolution,  rapid 
data  acquisition,  relatively  easy  interpretation  of  data,  capability  for 
at  least  semiquantitative  analysis,  wide  elemental  detection  range  (all 
elements  except  H and  He)  and  high  sensitivity  for  low  atomic  number 
elements  commonly  used  in  glass  constituents. 

Auger  Electron  Spectroscopy  is  generally  considered  to  be  a "non- 
destructive" technique.  However,  there  are  many  instances  where  the 
incident  electron  beam  can  alter  the  surface  composition  of  the  solid 
being  analyzed  [7-67].  Compositional  changes  in  glass  during  electron 
irradiation  for  AES  have  been  reported  by  several  investigators  [11,35, 
36,39,44,45,48].  Changes  are  particularly  severe  in  glasses  containing 
alkali  ions.  An  increase  or  decrease  of  intensity,  a change  in  energy 
and/or  a change  in  Auger  peak  shape  may  be  associated  with  electron  bom- 
bardment time.  Such  effects  must  be  at  least  minimized  in  order  to 
perform  a meaningful  analysis. 

The  purpose  of  this  research  is  to  obtain  a better  understanding  of 
the  effects  of  electron  beam  irradiation  upon  the  surface  of  glasses. 

The  topics  are  rather  specific  to  the  problem  associated  with  sodium  ion 
diffusion  and  desorption  in  glasses  induced  by  electron  irradiation  for 
Auger  electron  analysis,  but  the  conclusions  will  be  shown  to  be  generally 
applicable  to  other  alkali  ions. 
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Background  and  Qualitative  Assessment 

Interaction  of  Electrons  With  Glass 

When  a glass  is  exposed  to  an  electron  beam,  numerous  changes  in 
the  physical  and  chemical  properties  of  the  glass  surface  can  take  place. 
For  example,  the  glass  surface  composition  can  be  altered,  an  optical 
absorption  band  may  be  induced,  the  density  of  the  glass  surface  can 
change,  paramagnetic  defect  centers  can  be  formed,  microscopic  phase 
separation  can  be  induced  and  the  physical  and  chemical  strength  of  the 
glass  surface  can  be  altered  [68,69].  These  phenomena  all  result  from 
the  interaction  between  electrons  and  the  glass  surface  and  will  be 
referred  to  collectively  as  "electron  beam  effects." 

Electron  beam  effects  are  consequences  of  interaction  of  electrons 
with  the  glass.  Impinging  electrons  with  energies  less  than  20  keV  are 
scattered  elastically  or  inelastically  [70]  and  can  cause  changes  in  the 
electronic  state  of  atoms  near  the  surface.  The  electrons  move  through 
the  glass  matrix  and  will  either  create  a dangling  bond  in  the  glass 
structure  [68,69],  be  trapped  by  pre-existing  flaws  to  form  defect  centers 
in  the  glass  structure  [68],  recombine  with  the  positively-charged  hole 
[68],  or  produce  a secondary  electron  cascade  by  knock-on  collisions  with 
the  bound  electrons  [68,69].  Secondary  electrons,  created  by  primary 
electrons  through  the  matrix,  ionize  additional  bound  electrons  by 
Coulombic  interactions  and  lose  energy  for  each  ionization.  Finally, 
the  electrons  will  not  have  enough  energy  to  ionize  the  atom,  therefore 
they  become  trapped  at  defects  or  recombinewith  positive  holes  in  the 
matrix.  These  electronic  excitations  result  in  damage  near  the  surface. 
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Electron  Beam  Effects 

Although  the  early  interest  in  the  radiation  effect  on  glass  ori- 
ginated from  the  observation  by  Roentgen  that  the  glass  door  knob  on  the 
inside  of  his  x-ray  room  turned  blue,  a realization  that  electron  beams 
affect  glass  is  relatively  new.  Such  effects  were  first  observed  on 
glass  envelopes,  which  limited  the  lifetimes  of  cathode  ray  tubes  [49]. 

Todd  and  Lineweaver  [50]  found  that  the  electron  bombardment  of  glasses 
containing  alkali  elements  led  to  oxygen  outgassing  and  browning  of  the 
glass.  These  phenomena  have  usually  been  explained  in  terms  of  a mechan- 
ism proposed  by  Lineweaver  [51].  According  to  his  mechanism,  the  high 
energy  electrons  traverse  a thickness  equivalent  to  their  range  and  come 
to  rest,  producing  a net  negative  charge  in  the  glass.  This  net  electrical 
charge  creates  an  electric  field  in  the  glass.  Under  the  influence  of 
this  electric  field,  the  positive  alkali  ions  present  in  the  glass  diffuse 
away  from  the  surface  and  the  oxygen  ions  move  toward  the  surface  and  are 
desorbed  into  the  vacuum. 

Several  Electron  Microprobe  (EMP)  studies  on  alkali-containing 
glasses  [52-61]  apparently  supported  the  Lineweaver  mechanism.  The  decay 
of  the  x-ray  signal  from  the  alkali  atom  could  only  mean  that  alkali  atoms 
were  leaving  the  volume  in  which  the  incident  electrons  had  sufficient 
energy  to  excite  x-ray  emission.  Although  the  creation  of  an  internal 
electric  field  under  electron  bombardment  is  likely,  no  measurement  has 
been  reported  on  the  magnitude  of  the  electric  field. 
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Compositional  Change  in  Auger  Electron  Spectroscopy 

In  the  1970' s,  Auger  Electron  Spectroscopy  (AES)  became  popular 
and  was  applied  to  glass  surface  analysis.  Goldstein  and  Carlson  [62] 
measured  Auger  spectra  for  an  alkali-strontium-silica  glass,  while  Rynd 
and  Rastogi  [63,64]  used  AES  to  characterize  glass  fiber  surfaces. 

Chappell  and  Stoddart  [35]  used  AES  to  investigate  the  effects  of  the 
float  glass  process  upon  the  surface  composition.  Pantano,  Clark  and 
Hench  [65]  used  AES  to  examine  the  sequence  of  chemical  reactions  at  a 
soda-lime-silica  glass  surface  during  aqueous  attack.  A common  observa- 
tion in  the  above  studies  was  the  absence  of  alkali  at  the  surface  or 
the  rapid  decrease  of  alkali  Auger  signal  in  a fashion  similar  to  the 
EMP  study.  However,  the  problem  in  the  case  of  AES  is  much  more  severe, 
since  the  Auger  escape  depth  is  ~ 10  A [66]  in  comparison  to  ~ 1 pm  in  the 
EMP  [67].  Therefore,  any  change  within  a few  atomic  layers  causes  a 
noticeable  variation  in  the  Auger  signal  intensity. 

Figure  1(a)  and  (b)  compare  the  Auger  spectra  of  a freshly  created 
surface  and  a surface  with  extensive  electron  irradiation  prior  to  Auger 
measurement  of  soda-disilicate  glass.  Obvious  differences  are  seen  in 
the  intensity  of  the  sodium  Auger  peaks.  Comparing  the  spectrum  in 
Figure  1(b)  with  an  Auger  spectra  from  fused  silica  shown  in  Figure  1(c), 
little  difference  is  noticed  in  terms  of  peak  height  ratio  of  Si/0.  This 
indicates  that  electron  irradiation  changes  the  surface  composition  of 
the  soda-disilicate  glass  to  a composition  similar  to  fused  silica. 

This  similarity  is  further  substantiated  by  infrared  reflection 
spectra,  as  shown  in  Figure  2.  Curves  (1),  (2)  and  (3)  are  the  spectra 
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ELECTRON  ENERGY  IN  eV 

Figure  1.  Comparison  of  AES  spectra  for  (a)  a freshly  created  surface 
and  (b)  a surface  with  extensive  electron  irradiation  prior 
to  the  measurement  of  soda-disilicate  glasses,  and  (c)  a 
fused  silica.  A 3 keV,  5 mA/cm^  electron  beam  was  used. 
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Figure  2.  Comparison  of  IRRS  spectra  for  soda-disilicate  glass  surface  before  (curve  (1)) 

and  after  (curve  (2))  electron  bombardment,  and  a fused  silica  surface.  A 25  keV 
electron  beam  was  used  for  bombardment  in  a scanning  electron  microscope  (JEOL-35c) 
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from  the  soda-disilicate  glass  surface  before  and  after  electron  beam 
irradiation  and  from  the  fused  silica,  respectively.  The  spectrum  for 
the  electron  irradiated  surface  shows  a shift  in  the  Si-0  vibrational 
mode  toward  the  value  observed  for  the  fused  silica  (1100  cm  ■*■)  [71]. 

The  Si-0  vibrational  mode  in  curve  (2)  also  assumes  an  overall  shape 
similar  to  that  observed  from  fused  silica.  A NBO  stretching  peak  at 
960  cm  ^ [71]  observed  in  curve  (2)  is  very  weak  as  compared  to  curve 
(1) . This  indicates  that  the  sodium  is  depleted  at  least  within  the 
infrared  information  depth  (~  .5  pm  [72]). 

The  alkali  Auger  signal  loss  can  be  caused  by  the  disappearance  of 
alkali  at  the  surface  due  to  either  its  diffusion  into  the  bulk  or  its 
escape  into  the  vacuum.  Surface  diffusion  can  also  cause  the  composi- 
tional change,  particularly  in  Auger  microanalysis  when  an  extremely 
small  electron  beam  is  used  [ 73] . In  a typical  analysis  of  a glass 
surface,  either  an  electron  beam  scanning  a relatively  large  area  or  a 
de-focused  electron  beam  can  be  used  to  reduce  the  electron  beam  current 
density.  When  this  is  done,  surface  diffusion  can  be  ruled  out  because 
a measurable  signal  loss  is  greater  than  the  radius  of  the  electron  beam. 
Since  a measurable  signal  loss  still  occurs,  both  diffusion  into  the  bulk 
and  escape  from  the  surface  has  to  be  studied  more  carefully  in  order  to 
understand  the  beam  effects. 

Electron  Beam  Induced  Alkali  Diffusion 

The  surface  of  a glass  is  normally  electrically  isolated  from  ground 
but  has  an  incident  electron  current  entering  and  a secondary  electron 
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current  leaving  the  sample.  The  electrical  state  of  the  surface  is 
therefore  largely  governed  by  the  ratio  of  outgoing  to  incoming  current, 
i.e.,  by  the  secondary  electron  emission  yield,  5.  The  dependence  of 
the  yield,  6,  on  primary  beam  energy,  Ep,  is  shown  in  Figure  3 [70]. 

When  the  yield  is  less  than  unity,  the  surface  is  negatively  charged  due 
to  the  accumulation  of  negative  charges  on  the  surface.  Under  these 
conditions,  meaningful  Auger  spectra  are  often  difficult  to  obtain  [74]. 
The  surface  potential  may  reach  a high  negative  value  close  to  the  break- 
down voltage  and  dielectric  breakdown  may  occur  [75]  and  disturb  the 
Auger  spectra.  A positive  surface  potential  develops  when  the  yield  is 
greater  than  unity.  This  occurs  in  the  energy  range  between  E ^ and  E^ 
where  most  AES  data  are  collected.  In  this  case,  low  energy  secondary 
electrons  can  return  to  the  surface  and  prevent  excessive  surface  poten- 
tials. The  higher  6 is,  the  higher  the  positive  charge  should  be.  In 
the  case  of  SiC^,  the  maximum  of  6 is  reached  near  400  eV  [76]. 

The  influence  of  the  angle  of  incidence  of  the  electron  beam  on  the 
surface  charge  is  also  critical.  When  this  angle  is  increased,  the  yield 
increases  and  the  positive  surface  charge  is  expected  to  increase.  This 
is  demonstrated  by  data  from  Jassen  et  al.  [77]  on  MgO.  The  oxygen  Auger 
peak  shifted  10  eV  lower  in  energy  as  the  angle  of  incidence  was  increased 
from  0°  to  80°  with  respect  to  the  plane  of  the  surface.  The  electron 
beam  energy  and  angle  of  incidence  are  thus  critical  to  prevent  charging 
[78].  Upon  development  of  a stable  positive  surface  potential,  the  Auger 
spectrum  can  be  obtained.  However,  the  positive  surface  potential  shifts 
the  Auger  peaks  anywhere  from  2 to  10  eV  lower  in  energy.  In  this  case, 
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Figure  3.  Variation  of  the  secondary  electron  emission  yield  as  a 
function  of  the  primary  beam  energy. 
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an  electric  field  is  produced  normal  to  the  surface  as  a result  of  the 
surface  positive  potential  and  the  diffused  negative  layer  at  a depth 
equal  to  the  electron  range.  This  condition  is  shown  schematically  in 
Figure  4.  The  electric  field  strength  can  be  estimated  by  the  energy 
shift  of  the  Auger  peak  divided  by  the  electron  range,  i.e.,  10  eV/1000  A 
= 10  V/cm  [79].  From  the  electrical  conduction  theory  [ 80],  there  is 
little  doubt  that  such  a high  field  can  lead  to  diffusion  of  alkali  ions 
in  the  glass.  The  formation  of  a sodium  deficient  layer  during  10  keV 
electron  bombardment  of  a soda-lime-silica  glass  was  demonstrated  using 
SIMS  by  Gossink  and  Lommen  [81].  Their  data  indicate  that  the  near 
surface  layer  is  deficient  in  sodium  and  is  bounded  by  a steep  rise  of 
the  sodium  signal  at  a larger  depth.  The  thickness  of  this  layer  in- 
creased with  electron  irradiation  time  until  it  became  close  to  the 
expected  maximum  range  of  electron  penetration.  These  data,  along  with 
the  results  of  the  infrared  analysis  presented  in  Figure  2,  are  direct 
evidence  that  field-enhanced  diffusion  of  alkali  ions  does  occur  into 
the  bulk  of  glasses.  Therefore,  these  data  support  the  Lineweaver 
mechanism,  at  least  in  a qualitative  sense. 

As  reported  above,  a number  of  workers  [10,11,33,36,39,44,45,48] 
have  reported  field-enhanced  diffusion  of  mobile  ions  during  AES  analysis. 

For  example,  Chou  et  al.  [33]  observed  the  diffusion  of  Cl  to  the  surface 
of  thin  oxide  on  Si  during  AES  analysis.  Dawson  et  al.  [44]  observed 
time  dependent  change  in  Na,  Ca,  K,  C and  0 Auger  signals  from  medieval 
glasses.  They  attributed  the  change  in  Na,  K and  Ca  to  field-enhanced 
migration  of  cations.  Although  the  simple  mechanism  postulated  by  Lineweaver 
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Figure  4.  Schematic  representation  of  the  electric  field  produced 

normal  to  the  surface  as  a result  of  the  surface  positive 
potential  and  the  diffused  negative  layer  at  a depth  equal 
to  the  electron  range. 
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is  attractive,  neither  an  experimental  determination  of  the  field  strength 
nor  a mathematical  model  for  field-enhanced  diffusion  has  been  reported. 

Pantano  et  al.  [39]  developed  an  AES  technique  to  improve  measurement 
of  alkali  concentrations  at  glass  surfaces.  They  postulated  that  the 
alkali  disappeared  from  the  surface  by  diffusion.  The  rate  of  disappear- 
ance is  therefore  proportional  to  the  diffusive  flux  which  is  equal  to 
the  alkali  mobility  multiplied  by  the  driving  force  for  diffusion.  They 
postulated,  therefore,  that  the  high  mobility  of  sodium  in  glass  (which 
allowed  diffusion  away  from  the  region  of  investigation)  could  be  lowered 
by  reducing  the  temperature  of  glass.  Because  of  lower  mobility,  lower 
signal  decay  rate  would  be  observed.  Their  data,  reproduced  in  Figure  5, 
are  consistent  with  their  postulate.  An  additional  significant  result 
is  that  a characteristic  "incubation  time"  is  observed  prior  to  the  decay 
at  low  temperatures.  They  interpreted  the  incubation  time  on  the  basis 
of  a time  delay  before  the  local  temperature  (caused  by  electron  beam 
heating)  increased  to  a critical  temperature,  T , where  the  alkali  mobility 
was  sufficiently  high  for  field-enhanced  diffusion.  Although  there  is  no 
doubt  that  the  incident  electron  beam  causes  heating,  a number  of  experi- 
mental and  theoretical  studies  of  beam  heating  [ 82-  86  ] show  an  instan- 
taneous rise  in  temperature.  Therefore,  it  seems  unlikely  that  it  would 
take  so  long  for  the  temperature  to  rise  to  the  critical  temperature. 

For  example,  Mullendore  et  al.  [ 87  ] reported  that  the  AES  electron  beam 
caused  the  temperature  of  loose  titanium  powder  to  rise  to  the  critical 
temperature  for  oxygen  diffusion  (~  500°  C)  in  less  than  one  second.  Thus, 
a reasonable  explanation  for  the  incubation  time  still  remained  to  be 
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Figure  5.  Normalized  Auger  sodium  signal  decay  from  a soda-lime-silica 

glass  at  both  room  temperature  and  liquid  nitrogen  temperature 
using  a beam  diameter  of  500  pm  at  3 keV  for  various  levels  of 
beam  current  (from  [39]). 
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developed.  This  will  be  accomplished  in  the  present  work.  Nevertheless, 
the  incubation  time  effects  allow  the  use  of  AES  as  an  analytical  tech- 
nique for  structural  analysis  of  glass  and  ceramics,  as  demonstrated  in 
this  work. 

Electron  Beam  Induced  Alkali  Volatilization  and  Desorption 

In  considering  the  Auger  signal  decay  of  alkali  in  glass  observed 
with  AES,  one  must  also  consider  alkali  volatilization  and  desorption. 

This,  in  a sense,  is  the  opposite  of  diffusion  into  the  bulk.  Insulating 
materials  are  expected  to  be  most  susceptible  to  beam  heating.  Likewise, 
local  beam  heating  is  especially  prevalent  in  Auger  microanalysis  when 
high  beam  energies  and  small  beam  sizes  are  used.  The  examination  of 
small  or  thin,  unsupported  glass  specimens  presents  the  most  severe  beam 
heating  problems  because  of  the  difficulty  associated  with  removing  the 
energy  deposited  by  the  electron  beam.  Thus,  the  analysis  of  thin  glass 
membranes  was  unsuccessful  due  to  the  melting  of  holes,  as  shown  in  Figure 
6(a).  Another  example  is  given  in  Figure  6(b),  where  a small  glass  fiber, 
which  was  suspended  by  one  end,  was  observed  to  soften  and  collapse  when 
irradiated  by  an  electron  beam. 

Melting  of  glass  membranes  and  bending  of  glass  fibers  are  obviously 
extreme  cases.  However,  these  examples  indicate  that  the  electron  beam 
can  raise  the  temperature  of  the  sample  higher  than  the  liquidus  tempera- 
ture of  the  membrane  and  higher  than  the  softening  point  of  the  fibers. 

From  these  results,  there  is  little  doubt  that  electron  irradiation  can 
induce  alkali  volatilization  when  an  electron  beam  with  extremely  high 
current  density  is  used,  or  when  a specimen  is  isolated  from  a good  heat  sink. 
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Figure  6.  a)  Optical  microscope  photograph  of  a soda-lime-silica 
glass  membrane  showing  a hole  by  melting  as  a result  of 
electron  beam  heating,  b)  A photograph  showing  collapsed 
soda-lime-silica  glass  fiber  as  a result  of  electron  beam 
heating.  A 3 keV,  100  mA/cm^  electron  beam  was  used. 
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Another  consequence  of  electron  bombardment  is  electron  stimulated 
desorption  (ESD) . It  is  known  that  [70]  an  electron  beam  can  cause 
desorption  of  atoms  and  molecules  from  a surface  as  (1)  neutral  molecules, 
atoms  or  fragments,  (2)  excited  neutrals  or  (3)  positive  or  negative  ions. 
ESD  is  often  discussed  only  in  terms  of  adsorbates [88-92]  and  research 
on  desorption  from  a solid  is  limited  [93,94]. 

The  usual  picture  of  Electron  Stimulated  Desorption  according  to 
Menzel  and  Gomer  [95]  and  Redhead  [96  ] (MGR  model)  is  that  the  incident 
beam  excites  a bonding  electron  to  an  anti-bonding  state,  and  as  a result, 
the  effective  potential  between  surface  ions  and  the  solid  becomes  repul- 
sive and  the  ion  desorbs.  This  model  explains  many  of  the  distinct 
characteristics  of  ESD  observed  for  adsorbates  on  metal  surfaces.  The 
rate  of  electron  stimulated  desorption  from  adsorbates  is  given  by  [70]: 


de  -£e! 

dt  ' e 9 


(1) 


where  0 is  the  surface  coverage  in  the  binding  state  which  can  be  desorbed 
by  electron  impact,  a is  the  total  desorption  cross  section,  ip  is  the 
electron  current  density  and  e is  the  electronic  charge.  Therefore,  the 
coverage  exponentially  decreases  with  electron  bombardment  time  as  follows 


0 = exp(-t/x) 

where  (2) 

t = e/j^o  . 


In  glasses,  the  electropositive  alkali  atoms  are  bound  ionically 
to  non-bridging  oxygen  [ 80  ] . Applying  the  MGR  model  to  alkali  desorption 
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from  a glass  surface,  the  oxygen-sodium  interaction  potential  has  to 
become  repulsive  and  it  is  necessary  for  the  oxygen,  which  is  origi- 
nally ionized  negatively,  to  become  positively  charged.  This  means  that 
at  least  two  or  three  electrons  have  to  be  transferred  since  the  oxygen 
is  typically  bonded  with  a charge  between  -1  and  -2.  Therefore,  the  MGR 
model  cannot  explain  how  ESD  can  occur  in  glasses. 

Knotek  and  Feibelman  [ 97  ] have  recently  developed  a new  mechanism 
for  ESD,  namely,  an  interatomic  Auger  de-excitation  process.  This  mech- 
anism can  explain  the  large  charge  transfer  involved  in  electron  stimulated 
desorption.  The  desorption  mechanism  is  initiated  by  production  of  a core 
hole  on  the  cation  or  the  anion  followed  by  an  interatomic  or  intraatomic 
Auger  process  and  desorption  of  0+  via  a Coulombic  repulsion.  The  loss 
of  oxygen  leads  to  a reduction  of  surface  oxide.  The  electronic  structure 
of  glass  is  similar  to  that  of  maximal  valency  oxide;  therefore,  alkali 
desorption  along  with  oxygen  desorption  is  expected.  However,  there  has 
been  no  report  of  alkali  desorption  from  glasses  until  the  present  study. 

Research  Plan 

The  above  review  shows  that  a fundamental  and  complete  understanding 
of  mechanisms  of  alkali  Auger  signal  loss  has  not  yet  been  achieved.  This 
includes  the  theory  of  alkali  Auger  signal  decay.  Therefore,  the  present 
study  addresses  the  following: 

1.  The  parameters  describing  the  rate  of  alkali  Auger  signal  decay 
will  be  clarified  and  their  relative  magnitudes  will  be  discussed. 

2.  Consistent  models  of  the  process  will  be  developed  mathematically. 
The  experimental  results  will  be  compared  to  the  model. 
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In  Chapter  II,  changes  in  the  alkali  Auger  signal  as  a result  of 
electron  beam  effects  were  measured  for  thin  glass  films  deposited  on 
metal.  With  these  samples,  the  effects  on  driving  force,  mobility  and 
electron  stimulated  desorption  were  studied.  These  measurements  were 
coupled  with  theoretical  modeling  of  alkali  diffusion  with  emphasis  on 
the  field-enhanced  diffusion. 

Although  electron  stimulated  desorption  has  been  observed  for  a 
number  of  gas  adsorbates,  very  little  is  known  about  desorption  of  metal 
ions  from  oxide  surfaces.  In  particular,  there  is  a total  lack  of  data 
for  alkali  desorption  from  glass  surfaces.  Hence,  Chapter  III  is  devoted 
to  the  investigation  of  ESD  with  alkali  ions  during  analysis  with  AES. 
Crucial  experiments  were  designed  to  study  ESD.  A mathematical  model 
was  developed  to  explain  the  ESD  kinetics  and  the  results  of  the  experi- 
ments are  discussed  quantitatively.  The  mechanism  for  ESD  of  sodium 
from  soda-silica  glass  was  also  investigated  and  a model  for  sodium 
desorption  was  proposed. 

On  the  basis  of  the  work  presented  in  Chapters  II  and  III,  a general 
mathematical  description  of  the  Auger  signal  decay  for  mobile  ions  is 
presented  in  Chapter  IV.  The  theory  is  applied  to  low  temperature  measure- 
ment, and  the  origin  of  the  incubation  time  is  discussed.  Present  experi- 
mental results  are  compared  with  this  model. 

Finally,  in  Chapter  V,  it  is  shown  that  the  electron  beam  effects 
can  be  used  to  provide  a basis  for  Auger  electron  spectroscopy  of  micro- 
structural  as  well  as  compositional  analysis.  In  the  present  study, 
binary  lithia-disilicate  glass  was  chosen,  where  controlled  nucleation- 
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crystallization  of  a glass  provides  a microstructure  of  uniform  small 
crystals  surrounded  by  a glass  matrix.  The  composition  of  the  crystal 
is  the  same  as  that  of  bulk  glass  for  this  stoichiometric  material. 

Due  to  the  residual  phase,  compositional  variation  via  any  of  the  mech- 
anisms described  may  also  occur  in  glass-ceramics  during  electron  bom- 
bardment. The  electron  beam  effects  on  the  composite  material  are 
discussed. 


CHAPTER  II 

DIFFUSION  OF  SODIUM  DURING  ELECTRON  BOMBARDMENT 


Bombardment  of  the  glass  by  an  electron  beam  will  break  bonds  to 
form  ion-electron  pairs  [98]  and  will  raise  the  local  temperature  of 
the  material  [ 82-86  ].  Both  effects  cause  an  increase  in  the  alkali 
ion  mobility,  which  increases  the  rate  of  alkali  diffusion.  While  the 
mobility  of  the  alkali  ion  is  important,  a driving  force  is  necessary 
for  diffusion,  and  this  is  provided  by  an  electric  field.  If  the  charge 
build-up  inside  the  glass  upon  electron  bombardment  is  varied  by  some 
technique,  alkali  ion  diffusion  will  change.  To  change  the  space  charge, 
glass  thin  films  were  deposited  on  a conducting  material  so  that  energetic 
electrons  penetrated  through  the  film  and  were  stopped  by  the  conducting 
substrate.  Therefore,  it  was  postulated  that  the  resultant  electric 
field  would  be  changed  when  the  film  thickness  was  varied.  Electron 
beam  heating  would  also  be  minimized  since  good  contact  to  a thermal 
sink  was  provided.  With  these  samples,  the  region  of  the  glass  over 
which  forces  may  act  and  diffusion  will  occur  is  well  defined.  There- 
fore, a mathematical  model  can  be  developed  and  compared  with  the  results 
to  determine  important  parameters  such  as  the  electric  field. 

Experimental  Procedure 

O O 

The  measurements  reported  in  this  Chapter  were  from  1000  A or  2000  A 
soda-silica  glass  films  deposited  on  optically  polished  stainless  steel. 
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High  purity  SiC^  and  reagent  grade  ^2^^  were  simultaneously  evaporated 
from  an  electron  beam-heated  crucible  and  a resistively-heated  crucible, 
respectively.  During  deposition,  vacuum  was  maintained  at  5 x 10  ^ Torr 
and  the  substrate  was  heated  to  300°  C.  Film  thicknesses  were  controlled 
with  a quartz  crystal  monitor. 

For  comparison  with  thin  films,  bulk  glass  of  Na20-2Si02  was  prepared 
by  melting  mixtures  of  reagent  grade  ^2^^  and  high  purity  Si02  at  1400  ° C 
for  24  hours  in  a Pt  crucible.  After  casting  into  a graphite  mold,  it  was 
annealed  at  250°  C for  four  hours. 

The  Auger  measurements  were  made  in  a Physical  Electronics  Thin  Film 

Analyzer^  or  a Varian  Auger  Spectrometer"^  in  a residual  vacuum  pressure 

“8 

of  less  than  1 x 10  Torr.  All  the  samples  were  irradiated  by  an  elec- 
tron beam  at  an  angle  of  30  ° (with  respect  to  the  plane  of  the  surface) 
to  prevent  charging. 

In  the  former  instrument  (PHI) , electron  beams  with  a diameter  of 
600  ym  (defocused  beam)  and  energies  up  to  3 keV  were  used.  The  samples 
were  mounted  on  a carousel  which  was  cooled  with  liquid  nitrogen.  The 
specimen  temperature  was  monitored  with  a Chromel-Alumel  thermocouple^"^ 
sandwiched  between  thin  glass  slides.  The  Varian  instrument  was  capable 
of  producing  an  electron  beam  with  a diameter  of  10  ym  and  energies  up 
to  10  keV . A scanning  electron  beam  was  used  to  more  accurately  control 
the  electron  current  density  in  the  Varian  spectrometer.  The  total  beam 


t CMA  10-155,  Perkin  Elmer,  Physical  Electronics  Division,  Eden  Prairie,  MN 
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current  was  measured  by  a positively  biased  Faraday  cage  and  beam  current 
densities  were  calculated  based  upon  the  area  scanned.  Sputter  profiling 
was  accomplished  using  2 keV  Ar+  with  the  vacuum  chamber  backfilled  to 
2 x 10  ^ Torr. 

In  studying  the  change  in  sodium  Auger  peak  height  as  a function  of 
electron  beam  irradiation  time,  the  KLL-sodium  peak  at  990  eV  was  moni- 
tored by  scanning  the  energy  range  of  940  to  1020  eV.  A derivative 
signal (EdN/dE)  was  continuously  recorded  on  a strip-chart  recorder. 

Results 

O 

The  Auger  spectrum  from  1000  A thick  glass  film  deposited  on  a stain- 
less steel  substrate  is  shown  in  Figure  7.  The  spectrum  was  taken  using 
a 3 keV  electron  beam  with  a total  beam  current  of  5 pA  at  a liquid  nitro- 
gen temperature  (~  110 ° K) . A modulation  voltage  of  5 V peak-to-peak 
was  used.  The  KLL-sodium  peak  was  again  measured.  The  peak  height  had 
significantly  decreased,  as  shown  in  Figure  7.  No  elements  other  than 
Si,  C,  0 and  Na  were  detected.  Using  the  same  procedure,  several  spots 
were  analyzed  on  the  film.  It  was  found  that  the  relative  peak  heights 
for  Si,  0 and  Na  varied  by  less  than  5%  from  different  analysis  sites. 
Although  it  is  difficult  to  accurately  quantify  the  Auger  data,  the  glass 
film  composition  was  estimated  (using  bulk  glass  standards)  to  be  15  ± 5 
mole  % of  Na20  and  85  ± 5 mole  % of  Si02* 


Sputter  Ion  Gun  04-161,  Perkin  Elmer,  Physical  Electronics  Division, 
Eden  Prairie,  MN 
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ELECTRON  ENERGY  IN  eV 

Figure  7.  AES  spectrum  from  a 1000  A thick  glass  film  deposited  on 
a stainless  steel  substrate,  taken  using  a 3 keV  electron 
beam  with  a total  current  of  5 pA  at  liquid  nitrogen  tem- 
perature (110°  K) . 
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The  peak  heights  (normalized  to  the  peak-to-peak  height  at  zero 

time)  for  the  sodium  signal  are  plotted  as  a function  of  beam  irradia- 

-3  2 

tion  time  in  Figure  8.  A3  keV,  5 x 10  A/cm  electron  beam  was  used 

° ° 

to  excite  the  sodium  Auger  electrons  from  1000  A and  2000  A films  and 
from  bulk  glass.  The  measurements  were  made  at  room  temperature.  The 
rate  of  sodium  signal  decay  was  strongly  dependent  upon  the  film  thick- 
ness, with  significantly  slower  decay  observed  for  thin  films  as  compared 
to  the  bulk  glass.  The  times  required  for  the  signal  to  decay  to  50% 
of  its  initial  value  were  approximately  5.5,  2.8  and  0.4  minutes  for  the 
1000  A and  2000  A glass  films  and  the  bulk  glass,  respectively.  When  a 
beam  current  density  of  less  than  3 mA/cm  was  used,  the  sodium  signal 
first  increased  for  one  minute  and  then  decayed.  Dawson  et  al.  [44]  made 
a similar  observation  and  attributed  the  effects  to  removal  of  a carbon 
overlayer  which  inelastically  scattered  the  sodium  Auger  electrons.  In 
the  present  case,  the  carbon  signal  was  observed  initially  (Figure  7), 
but  it  decreased  rapidly  with  electron  bombardment.  The  effect  of  removal 
of  carbon  contamination  only  delays  the  time  to  start  sodium  Auger  signal 
decay,  but  the  changes  in  the  signal  decay  rate  versus  film  thickness  are 
the  same  as  those  observed  at  higher  beam  current  densities.  This  indi- 
cated that  carbon  contamination  on  the  surface  can  cause  significant  error 
in  evaluating  the  surface  composition.  Therefore,  minimal  surface  contam- 
ination^ was  desired  for  analysis. 

The  decay  rate  of  the  sodium  signal  was  also  dependent  upon  the 

O 

beam  energy.  The  time  dependences  of  the  Na  Auger  signal  from  a 2000  A 
thick  film  with  energies  ranging  from  1.5  to  10  keV  are  shown  in  Figure  9. 


Na  Auger  Signal  Decay  Curve 
Film  Thickness  Dependence 
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All  the  measurements  were  made  at  room  temperature.  In  these  measurements, 
a beam  current  density  of  1.4  x l(f3  A/cm2  was  used.  As  pointed  out 
above,  the  carbon  contamination  seriously  influenced  the  measurements  and 
the  increase  in  the  sodium  signals  for  the  first  few  minutes  resulted 
from  contamination  desorption.  Nevertheless,  the  dependence  of  the  sodi 
urn  Auger  signal  decay  upon  beam  energy  is  obvious.  For  the  1.5  keV  elec 
tron  beam,  the  sodium  decay  was  quite  rapid,  while  relatively  long  plateaus 
were  observed  for  3.0  and  7.5  keV  electrons.  Almost  no  decay  was  observed 
from  the  10  keV  beam.  These  differences  in  decay  rate  will  be  attributed 
to  changes  in  the  electric  field,  which  is  the  most  important  driving  force 
for  sodium  diffusion. 

The  results  obtained  for  2000  l glass  film  with  3 keV  electrons  when 

—3  2 2 

the  beam  current  density  was  varied  from  1.4  x 10  A/ cm  to  1.25  A/ cm  , 

are  shown  in  Figure  10.  As  the  current  density  was  increased,  the  decay 
rates  increased  drastically.  According  to  beam  heating  models  [82-86], 
the  steady  state  temperature  rise  is  directly  proportional  to  the  beam 
power,  which  is  the  product  of  energy  and  current.  At  a constant 
electron  energy,  the  temperature  and  therefore  the  mobility  of  sodium 
becomes  higher  as  the  current  density  is  increased.  Consequently,  faster 
diffusion  rate  should  result.  This  feature  is  clearly  seen  in  Figure  10, 
however,  there  is  no  simple  relationship  (e.g.,  linear,  quadratic,  etc.) 
between  the  rate  of  decay  and  the  beam  current  density. 

The  sodium  mobility  in  the  glass  can  also  be  changed  by  varying  the 
substrate  temperature.  The  sodium  Auger  signal  decays  at  110°  K and  300°  K 
are  compared  in  Figure  11.  Lower  temperatures  result  in  a slower  decay 


rate  when  other  factors  are  constant. 


Na  Auger  Signal  Decay  Curve 

Beam  Current  Density  Dependence  e- beam  = 3.0  KeV 
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Figure  10.  Sodium  Auger  peak  height  as  a function  of  beam  impingment  time  from  a 2000  A 

thick  film  at  room  temperature.  A 3 keV  electron  beam  was  used  with  the  current 
densities  varied  between  1.4  mA/cm^  and  1,25  A/ cm2. 


Na  Auger  Signal  Decay  Curve 
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Figure  11.  Sodium  Auger  peak  height  as  a function  of  beam  impingement  time  from  a 

1000  A thick  film  for  a 3 keV  electron  beam  with  current  density  of  5 mA/cm^ 
at  liquid  nitrogen  (110°  K)  and  room  temperature  (300°  K) . 
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A Model  for  Field-Enhanced  Diffusion 

For  thin  films,  the  sodium  signal  decay  rates  were  strongly  dependent 
upon  the  electron  beam  parameters.  The  decay  rate  can  be  controlled  by 
changing  either  the  sodium  mobility  in  the  glass  or  the  electric  field 
for  diffusion.  For  thick  glass,  energetic  electrons  enter  and  dissipate 
their  energy  by  ionization  or  other  electronic  excitation  of  the  atoms. 

They  come  to  rest  at  some  depth,  R,  within  the  glass  and  produce  a net 
negative  charge.  Theoretical  range  expression  can  be  obtained  by  inte- 
gration of  an  energy  loss  expression  [ 99  ] , which  is: 

R = •rE=E0  dE/dx  dE' 

A simplified  range  expression  was  developed  by  Cosselett  and  Thomas  [100], 

R = K Eq,  where  K is  a constant  and  n = 1.2  to  1.7,  depending  on  the  defi- 
nition of  range.  For  low  energy  electrons,  the  range  expression  of  Coss- 
lett  and  Thomas  is  not  precise  because  of  high  inelastic  scattering. 

However,  it  is  certain  that  as  the  primary  energy  is  increased,  the  elec- 
trons penetrate  further  into  the  solid.  Typical  ranges  are  a few  thousand 
Angstroms,  for  3 keV  electrons  and  more  than  one  micron  for  10  keV  electrons 
[101],  The  sodium  Auger  signal  decay  rate  was  strongly  dependent  upon  the 
film  thickness  for  a fixed  energy,  as  shown  in  Figure  8.  This  was  not 
true  in  the  bulk  glass,  where  the  decay  rate  was  almost  independent  of 
energy  when  the  current  density  was  kept  constant.  Therefore,  the  data 
shown  in  Figure  8 can  be  explained  in  part  by  postulating  that  as  the 
thickness  of  the  film  is  decreased,  the  amount  of  charge  trapped  also 
decreased,  causing  a lower  driving  force  and  a lower  diffusion  rate. 


The  dependence  of  the  decay  rate  upon  the  electron  energy  shown 
in  Figure  9 is  interpreted  in  the  same  manner.  This  is  shown  schemati- 
cally in  Figure  12.  When  1.5  keV  electrons  enter  the  2000  A glass  film. 


most  of  the  electrons  are  trapped  within  the  film  since  the  electron 
range  for  1.5  keV  is  almost  the  same  as  the  glass  film  thickness.  A 
space  charge  is  created  inside  the  film  and  diffusion  results.  For 
higher  energy  beams  (e.g.,  7.5  or  10  keV) , the  electrons  penetrate 
through  the  film  and  come  to  rest  in  the  metallic  substrate.  Therefore, 
the  number  of  electrons  trapped  inside  the  glass  film  at  high  energies 
is  significantly  lower  than  at  low  energies.  Correspondingly,  diffusion 
rates  are  slower  for  high  energy  electron  bombardment. 

The  discussion  presented  above  suggests  the  following  model  for 
sodium  diffusion  during  an  AES  analysis.  When  a space  charge  results 
from  bombardment  with  energetic  electrons,  the  negative  charge  will  be 
concentrated  at  a depth  about  equal  to  the  electron  range.  A positively 
charged  layer  will  also  exist  on  the  surface  since  the  secondary  electron 
emission  yield,  6,  is  normally  greater  than  unity.  A capacitor  structure 
can  be  used  to  model  the  sodium  diffusion  in  the  thin  film.  The  model 
consists  of  the  glass  thin  film  sandwiched  by  two  electrodes,  as  is 
illustrated  in  Figure  13.  The  diffusion  equation  including  the  electric 
field,  F (volt/cm),  becomes 


9c(x,t) 
9 t 


9F(x,t)c(x,t) 
P 9x 


(3) 


where  p is  the  mobility  of  the  sodium  and  D is  the  diffusion  constant. 
The  coordinate  perpendicular  to  the  surface  is  z.  For  the  present  case, 
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Low  Energy  High  Energy 


Schematic  illustration  of  the  electron  beam  penetration  and 
the  charge  accumulation  in  glass  film  for  low  and  high  energy 
electrons. 


Figure  12. 
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Figure  13.  Schematic  illustration  of  the  model  used  to 
calculate  sodium  ion  diffusion. 
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it  was  assumed  that  F is  independent  of  the  time  and  was  taken  as  an 
average  field  strength,  F . Then  equation  (3)  becomes 

clV 

9c(x,t)  = 32c (x, t)  _ 3c(x,t) 

at  „ 2 avp  3x  V ’ 

3x 

From  the  Einstein  relation  [102] the  mobility  of  the  sodium  is  expressed 
by 


y = eD/kT 


(5) 


where  e,  k and  T are  the  electron  charge,  Boltzman  constant  and  tempera- 
ture, respectively. 

Assuming  that  the  sodium  concentration  is  uniform  before  the  analysis 
(t  < 0)  leads  to  the  initial  condition 


c(z,t)  = cq  for  t < 0,  0 ■g  z < i 


(6) 


where  Cq  is  the  initial  concentration  of  sodium  in  the  film  and  a is  the 
film  thickness.  We  assume  that  no  electron  stimulated  desorption  is 
occurring  at  the  surface  and  sodium  is  physically  blocked  by  the  stain- 
less steel  at  the  glass/substrate  interface.  As  a result,  there  are  no 
fluxes  of  sodium  ions,  J,  at  the  surface  and  interface,  i.e., 


J = D — - F yc  — 0 for  x = 0 and  x = l. . 
3x  avH 


(7) 


Using  these  boundary  conditions,  the  equation  can  be  solved  analytically. 
The  solution  is  (see  Appendix  A) 

C0Favy£  exP<Favx/D) 


c(x,t)  = 


D(exp(Fy£/D)-l) 


+ c016D2FavyTT«,3  exp(Favyx/2D) 


» exp(  (4D2m2Tr2+F2vy2)  t / 4D2  g,2)m(  1-exp  (~Favyl/2P)  ) 


X 2 

m 


2mirD  nurx. 
— cos  — ), 


,/n2  2 2 2 2 2.2 

(4D  m it  +Fc,„y  Z ) 

cL V 


. mux 

(F  y sm  — - — + 
av  Z 


(8) 
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The  variation  of  sodium  concentration  within  the  film  after  electron 
bombardment  can  be  calculated  from  this  equation. 

The  sodium  surface  concentration  (i.e.,  at  z = 0)  as  a function  of 
electron  bombardment  time  is 


c(0,t)  _ 


F y£ 
av 


cQ  D(exp(Favy£/D)-l) 


+ 16D2F  yir£3 
av^ 


X l 
m 


2nnrD  exp(  (4D2m27r3+Favy3)  t/4D  £ ^)m (1-exp  (Favy£/2D)) 

9 2 2 2 2 2 2 

£(4D  i i +FZ  y £ ) 
av 


(9) 


Therefore,  the  normalized  Auger  signal  can  be  approximated  by  equation  (9), 
since  the  Auger  escape  depth  of  sodium  is  very  small  [66]  relative  to  the 
film  thickness. 


Discussion 


A Model  Calculation  and  Estimation  of  Field  Strength 

A theoretical  model  has  been  presented  for  describing  sodium  Auger 
signal  decay  in  the  glass  thin  film.  The  model  considers  only  diffusion 
influenced  by  the  electric  field  which  was  assumed  to  be  uniform  through- 
out the  film.  In  a quantitative  comparison  of  this  model  with  the  exper- 
mental  data,  the  field  strength,  F , can  be  treated  as  an  adjustable 
parameter  such  that  calculated  curves  give  a best  fit  with  experimental 
data.  As  discussed  before,  local  heating  by  the  electron  beam  causes 
the  temperature  of  the  substrate  to  increase,  thereby  increasing  the 
mobility.  This  beam  heating  effect  is  implicitly  included  in  the  model 
through  the  diffusion  constant. 
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As  has  been  shown  in  a number  of  studies  of  beam  heating,  the  local 
temperature  rises  almost  instantaneously  [ 82-86  ] and  reaches  a steady 
state  value  in  the  thin  film.  The  local  temperature  has  been  measured 
by  a thin  thermocouple  on  a glass  microscope  slide  and  temperature  rises 

O O 

of  150  to  250  C were  obtained  [103]  when  a 3 keV  electron  beam  with  a 

2 

beam  current  density  of  1 mA/cm  was  used.  Therefore,  for  analysis  with 

the  carousel  at  room  temperature,  a diffusion  constant  of  ~ 3.5  x 10  ^ 

2 

cm  /sec  (corresponding  to  the  constant  at  200°  C [104])  is  taken  as  an 
appropriate  value.  Calculated  results  from  equation  (9)  are  plotted  in 
Figure  14,  where  best  fits  to  the  experimental  data  presented  in  Figure 
8 are  given.  The  average  field  strengths  determined  from  these  calculations 
were  2.2  x 10^,  3.9  x 10^  and  1.2  x 10^  V/cm  for  1000  A and  2000  A films 
and  bulk  glass,  respectively.  For  the  bulk  glass,  the  electron  range 

O 

was  calculated  (3000  A for  3.0  keV  electron  beam  [101])  and  used  as  the 
equivalent  film  thickness.  Therefore,  the  estimated  field  strength  in 
the  bulk  glass  may  vary  by  a factor  of  one  to  five,  depending  on  the 
range  determined.  However,  the  calculated  field  strength  does  show  a 
consistent  trend  with  film  thickness. 

The  surface  charging  can  be  estimated  from  shifts  in  energy  of  the 
Auger  peak  [79].  For  the  2000  A film,  peak  shifts  of  1 ± 0.5  eV  were 

measured  for  a 3.0  keV  electron  beam.  Assuming  this  charge  acts  over 

° 4 

the  2000  A film,  the  field  strength  is  estimated  to  be  5 x 10  V/cm. 

This  is  approximately  the  same  as  the  field  strength  calculated  from  the 

model.  The  energy  of  the  sodium  Auger  transition  from  bulk  glass  shifted 

from  3 to  6 eV,  again  giving  a field  strength  of  105  to  2 x 10^  V/cm, 
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assuming  an  electron  range  of  3000  A;  this  is  reasonably  consistent  with 
the  value  derived  by  curve  fitting  the  experimental  results,  This 
coincidence  suggests  that  our  assumption  of  a uniform  average  field  is 
good  in  this  case.  However,  the  field  in  a bulk  sample  may  be  signifi- 
cantly different.  For  example,  the  electron  range  in  a bulk  film  is 
only  an  average  value  and  the  resulting  charge  will  straggle  over  a large 
distance  about  this  value  (e.g.,  ± 500  A for  a 3 keV  electron  in  Si02  [105]). 
In  addition,  a non-uniform  distribution  of  sodium  ions,  electrons  and 
holes  may  exist.  Therefore,  the  diffusion  equation  has  to  be  solved 
under  non-uniform  electric  field  conditions  (equation  (3)).  This  cannot 
be  done  until  a detailed  charge  distribution  is  known. 

Local  heating  by  the  electron  beam  has  been  approximated  in  this 
case  based  upon  experience  [103].  This  is  only  a crude  estimate.  Beam 
heating  was  expected  to  vary  not  only  when  the  current  density  was  changed, 
but  when  the  geometry  was  changed.  The  sodium  Auger  signal  decays  at 
110°K  and  300°K  (Figure  11)  were  not  significantly  different.  This  is 
because  the  film  is  in  direct  contact  with  a good  thermally  conducting 
material  and  heat  is  dissipated  well,  even  at  room  temperature.  The  temp- 
erature, thereby  the  mobility,  may  not  be  as  high  as  those  estimated  (i.e., 
200°C).  Therefore,  the  field  strength  in  the  film  should  be  higher  in 
order  to  provide  a larger  driving  force  and  diffusion  rate.  However,  the 
calculation  does  provide  a lower  limit  on  possible  field  strength  and  it 
indicates  that  a very  high  field  exists  in  the  film. 
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Sodium  Depth  Profile 

In  the  above  discussion,  the  sodium  Auger  signal  decay  was  depen- 
dent upon  the  electron  beam  energy,  current  density  and  substrate  tem- 
perature. It  was  found  that  the  field  induced  by  the  electron  beam  was 
about  105  V/cm.  In  such  a high  electric  field,  there  is  no  doubt  that 
the  sodium  will  move  [80,106].  According  to  the  model,  the  sodium  ions 
are  driven  into  the  thin  glass  but  are  physically  blocked  at  the  inter- 
face between  the  film  and  the  stainless  steel.  Thus,  accumulation  of 
the  sodium  at  the  interface  is  expected.  This  is  shown  by  the  concentra- 
tion versus  depth  curves  calculated  from  equation  (7).  The  calculated 
sodium  concentrations  versus  depth  in  the  film  after  irradiation  by  an 
electron  beam  is  shown  in  Figure  15.  Sodium  is  depleted  near  the  surface 
and  accumulates  at  the  interface.  This  accumulation  is  enhanced  by  a 
longer  irradiation  time. 

To  test  this  predicted  accumulation,  the  sodium  concentration  versus 

depth  was  measured  by  Ar"^"  sputtering.  A depth  profile  after  irradiation 

of  the  surface  of  a 1000  A glass  film  is  shown  in  Figure  16.  A3  keV, 

1 x 10  A/cm  electron  beam  was  scanned  over  ~ .2  mm  to  obtain  homogeneous 

irradiation  of  the  film.  After  the  signal  decayed  to  25%  of  its  maximum 

-3  2 

peak  height,  the  sodium  was  sputter  profiled.  A 3 keV,  1 x 10  A/cm 
primary  electron  beam  was  used  during  sputtering,  but  only  sodium  and 
iron  Auger  signals  were  recorded.  The  beam  bombardment  time  was  no  more 
than  10  seconds  at  each  depth  profile  concentration  point.  As  shown  in 
Figure  16,  lower  concentrations  of  sodium  were ^observed  in  the  middle  of  the 
film,  and  a large  accumulation  was  again  detected  at  the  interfa'ce.  The  depth 
profile  from  an  unirradiated  area  of  the  1000  A film  is  also  shown  for  the 


Concentration  C/C 
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Figure  15.  Calculated  sodium  distribution  in  the  film  after  electron 

beam  impingement  for  time  t,  when  t,  = 0 < t„  < t„  < t,  < t,-. 

i-  4 14  5 


Auger  Peak  Height 
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Figure  16.  Sodium  concentration  depth  profile  in  the 

1000  A glass  film  from  an  unirradiated  area 
(shown  by  dotted  curve)  and  an  irradiated 
area  (shown  by  solid  curve) . Iron  signal 
profile  is  shown  by  broken  line  curve. 
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comparison.  Higher  accumulation  at  the  interface  for  the  irradiated 
area  indicates  that  some  of  the  sodium  migrates  into  the  solid  due  to 
the  surface  electron  bombardment. 

However,  the  depth  profile  of  the  unirradiated  area  departs  signi- 
ficantly from  the  expected  uniform  concentration  distribution.  This  is 
thought  to  result  from  sodium  diffusion  induced  by  ion  bombardment. 

McCaughan  et  al.  [107]  observed  that  Ar+  bombardment  caused  migration  of 
sodium  ions  in  Si02*  An  accumulation  of  sodium  at  the  Si02/Si  interface 
was  observed  and  the  sodium  depth  profile  was  similar  to  Figure  16.  Al- 
though the  interaction  of  ions  with  the  solid  is  much  more  complicated  than 
that  of  electrons,  they  explained  the  result  by  considering  an  ion  neutrali- 
zation process  [108]  at  the  surface.  Incoming  low  energy  Ar"*"  are  neutralized 
within  a few  Angstroms  of  the  Si02  surface.  This  leads  to  a build-up  of 
a positive  charge  on  the  surface,  which  causes  sodium  diffusion.  In  the 
present  instance,  the  sodium  diffuses  more  than  the  Auger  escape  depth, 
resulting  in  a lower  Auger  signal  intensity.  The  glass/stainless  steel 
interface  blocks  sodium  diffusion  and  an  accumulation  results. 

In  each  depth  profiling  step,  the  electron  beam  was  exposed  for  a 
short  time  (less  than  10  seconds)  to  observe  the  sodium  concentration. 

In  this  case  the  electric  field  is  instantaneously  established  and  the 
sodium  would  move  into  the  bulk  under  the  field.  However,  the  estimated 
decay  within  10  seconds  is  only  ~ 3%  at  1 x 10  A/cm  , in  addition,  the 
sodium  may  relax  back  towards  the  surface  after  the  electron  beam  was 
turned  off.  Thus,  the  10  seconds  of  electron  beam  exposure  is  not  expected 
to  cause  a severe  depletion  of  sodium  in  the  middle  of  the  film. 
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In  spite  of  the  ion  beam  effect,  the  depth  profiles  in  Figure  16 
still  have  a discrepancy.  If  only  electron  beam  induced  diffusion 
occurred  during  Auger  measurements,  the  area  under  the  depth  profile 
should  be  proportional  to  the  total  amount  of  the  sodium  in  the  film. 

By  comparing  the  area  under  the  two  profiles,  it  is  apparent  that  the 
amount  of  sodium  ions  from  the  irradiated  area  (curve  (1))  is  smaller 
than  that  from  the  unirradiated  area  (curve  (2)).  This  indicates  that 
some  of  the  sodium  ions  are  lost  during  the  surface  irradiation,  i.e., 
electron  stimulated  desorption  has  occurred. 

In  this  case,  the  effect  of  electron  stimulated  desorption  should 

be  taken  into  account  in  the  model.  However,  as  will  be  shown  in  Figure 

28  (Chapter  III) , the  electron  stimulated  desorption  becomes  important 

2 

only  when  the  beam  current  is  high  (e.g.,  >1  mA/cm  ).  Thus,  the  field 
strength  calculated  from  the  present  model  is  still  valid.  The  importance 
of  electron  stimulated  desorption  is  substantiated  by  the  data  in  Figure 
10.  Anomalously,  high  decay  rates  observed  at  high  beam  current  densities 
(e.g.,  2 x 10  A/cm  or  1.25  A/cm  ) are  attributed  to  electron  stimulated 
desorption.  Therefore,  we  could  not  correlate  the  Auger  signal  decay  rate 
with  the  beam  current  density.  Both  field-enhanced  diffusion  and  electron 
stimulated  desorption  are  occurring  simultaneously  to  cause  the  sodium 
signal  decay.  Electron  stimulated  desorption  will  be  discussed  in  the 
next  chapter. 


Summary 

1.  Change  in  the  sodium  Auger  signal  due  to  electron  bombardment  has  been 
studied  using  thin  films  of  soda— silica  glass  deposited  on  a stainless 
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steel  substrate.  With  these  samples,  the  effects  of  electric  field, 
mobility  and  temperature  have  been  demonstrated. 

2.  The  amounts  of  charge  trapping  and  local  heating  are  a strong  func- 
tion of  the  beam  parameters  for  thin  films.  For  example,  the  time 
required  for  the  sodium  Auger  signal  to  decay  to  50%  of  its  initial 
value  increased  as  the  beam  energy  was  increased,  or  as  the  current 
density  was  decreased. 

3.  The  rearrangement  of  sodium  due  to  charge  trapping  was  modeled  using 
a configuration  similar  to  a capacitor.  A diffusion  equation,  inclu- 
ding the  electric  field,  was  solved.  The  calculated  and  experimental 
data  agree  well  and  indicate  that  fields  of  ~ 10^  V • cm  ^ exist 
during  analysis. 

4.  The  model  was  also  used  to  calculate  the  distribution  of  sodium  with 
depth,  but  the  calculated  distribution  disagreed  with  the  experimental 
sputter  profiling  data.  This  probably  results  from  redistributions 
during  ion  bombardment  for  sputter  profiling.  Therefore,  care  has  to 
be  taken  in  interpreting  depth  profiles. 

5.  Electron  stimulated  desorption  of  sodium  occurs  during  AES  analysis. 


CHAPTER  III 

DESORPTION  OF  SODIUM  DURING  ELECTRON  BOMBARDMENT 

As  stated  in  the  previous  chapter,  electron  stimulated  desorption 
(ESD)  is  one  of  the  mechanisms  for  alkali  Auger  signal  loss  during  analysis 
with  AES.  Upon  electron  bombardment,  some  of  the  sodium  atoms  are  desorbed 
into  the  vacuum  from  the  surface  and  others  diffuse  into  the  bulk  at  the 
same  time.  Therefore,  the  Auger  signal  decay  is  more  complicated,  since 
both  diffusion  and  desorption  have  to  be  considered  simultaneously.  The 
electron  stimulated  desorption  of  sodium  was  studied  in  the  AES-SIMS 
system  with  experimental  conditions  similar  to  those  of  conventional  AES 
analysis.  A phenomenological  model  was  developed,  in  which  diffusion  and 
desorption  are  simultaneously  taken  into  account.  The  results  are  then 
compared  with  the  model  and  an  ESD  cross  section  for  sodium  was  evaluated. 

A mechanism  for  sodium  desorption  from  glass  was  also  determined  by  applying 
a recently  developed  ESD  model  where  Auger  de-excitation  of  a core  hole 
causes  desorption  [ 97,109]. 

Experimental  Procedure 

A Na20-2Si02  glass  was  prepared  in  a manner  similar  to  that  described 
in  Chapter  II.  Glass  fibers  with  a diameter  of  400  to  800  ym  were  drawn 
from  the  molten  glass  state.  The  fibers  were  then  held  with  a clip  attached 
to  a carousel.  In  order  to  create  a new  clean  surface,  the  fibers  were 
fractured  in  the  vacuum  chamber  at  a pressure  of  1 x 10-^  Torr  just  prior 
to  the  ESD  experiment . 
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The  experimental  system  incorporated  a UTI  quadrupole  mass  filter^ 
with  a 3M  prefilter,  AES  single  pass  CMA,  a 3M  minibeam  ion  gun 
and  a sample  manipulator.  The  secondary  ion  mass  spectrometry  (SIMS) 

and  Auger  electron  spectrometry  (AES)  can  be  used  simultaneously  without 
changing  the  sample  location.  Therefore,  in  the  present  ESD  experiment, 
fractured  surfaces  were  exposed  to  the  scanning  electron  beam  from  the 
coaxial  gun  in  the  CMA,  and  desorbed  ions  were  detected  by  the  quadrupole 
mass  filter. 

In  studying  ESD  of  sodium  from  glass  surfaces,  the  quadrupole  mass 
filter  was  tuned  to  m/e  = 23  (Na"*")  and  the  intensity  was  recorded  as  a 
function  of  electron  bombardment  time.  During  measurements,  the  ion  gauge 
was  off.  The  background  signal  at  m/e  = 23  (such  as  Na  ions  from  the  W 
filament  of  the  electron  and  ion  guns,  Na  ions  from  the  vacuum  chamber 
walls,  etc.)  was  eliminated  by  setting  the  signal  discriminator 
at  the  background  signal  level. 

_ O 

Primary  electron  energies  of  500  to  4000  eV  and  beam  current  of  10 

_ 6 

to  10  A were  used  during  most  experiments,  since  these  are  the  ranges 
used  in  typical  Auger  analyses  [1] . The  electron  beam  current  density 
was  determined  by  dividing  the  total  beam  current  by  the  area  scanned. 
Auger  measurements  were  made  simultaneously. 


t UTI  100C,  UTI,  Sunnyvale,  CA 

ft  3M  Model  610  SIMS  Prefilter,  3M  Company,  St.  Paul,  MN 

CMA  15-110,  Perkin  Elmer,  Physical  Electronics  Div. , Eden  Prairie,  MN 
tttt  Minibeam  Ion  Gun,  3M  Company,  St.  Paul,  MN 

+++++  PM-600XYZTRC , Huntington  Mechanical  Lab,  Inc.,  Mountain  View,  CA 

tit+tt  Model  9302,  Amplifier-Discriminator,  0RTEC,  Inc.,  Oak  Ridge,  TN 
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A threshold  energy  measurement  for  desorption  of  sodium  was  carried 
out  by  rearranging  the  electronics  cabling  of  the  Auger  electron  spectro- 
scopy unit  and  the  mass  spectroscopy  unit,  as  shown  in  Figure  17.  The 
negative  sweep  voltage  normally  used  to  drive  the  cylindrical  mirror  ana- 
lyzer was  used  as  a sweepable  accelerating  voltage  for  the  electron  beam. 

In  this  case,  the  energy  of  primary  electrons  can  be  directly  read  from 

the  digital  panel  meter  housed  in  the  Auger  controller.^  The  beam  current 

”8 

was  maintained  in  the  10  A range  by  adjusting  the  emission  current  as 
the  beam  voltage  was  changed. 

Results 

Electron  stimulated  desorption  spectra  from  glass  surfaces  are  shown 

in  Figure  18(a)  and  (b).  The  sample  was  exposed  to  a pressure  of  (a) 

-9 

10  minutes  at  1 x 10  Torr  , and  (b)  90  minutes  after  an  ^n  sztu  fracture. 
The  quadrupole  mass  filter  was  scanned  from  m/e  = 1 to  50.  Major  peaks 
observed  in  the  spectra  are  1,  16,  17,  19,  23,  28  and  44  in  mass  to  charge 
ratio  (m/e),  which  are  identified  as  H+,  0+,  0H+,  F+,  Na+,  C0+  and  CC^"1”, 
respectively.  It  should  be  noted  that  all  peak  intensities  vary  with 
electron  beam  bombardment  time  and  the  relative  intensities  are  strongly 
dependent  upon  the  time  of  exposure  to  the  electron  beam.  Therefore,  the 
spectra  shown  in  Figure  18  are  only  representative  of  typical  data  but  do 
indicate  which  ions  are  representative  of  the  original  fracture  surface. 


t 11-500A  Auger  System  Control,  Perkin  Elmer,  Physical  Electronics  Division, 
Eden  Prairie,  MN 
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Figure  17.  Experimental  configuration  for  threshold  energy  measurement  of  ESD. 
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Figure  18.  ESD  spectra  from  a soda-silica  glass  surface 
exposed  to  a pressure  of  1 x 10“9  Torr  for  10 
minutes.  (a)  For  90  minutes,  (b)  after  an  in 
situ  fracture.  A 1 keV,  1 mA/cm^  electron  beam 
was  used. 
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The  relatively  high  peaks  at  m/e  =1,  17,  19  and  28  in  Figure  18(b) 
indicate  that  the  fracture  surface  is  partially  covered  with  H^O  and 
CO.  Partial  pressure  analysis  with  the  mass  filter  indicated  that  I^,  1^0 

_9 

and  CO  were  the  major  residual  gases  at  1 x 10  Torr  but  no  F compounds 
were  observed.  These  peaks  are  considerably  smaller  in  the  spectrum  taken 
10  minutes  after  the  fracture  (Figure  18(a))  and  peaks  for  m/e  = 1,  16 
and  23  are  dominant. 

This  observation  substantiates  that  a clean  surface  is  conserved  for 

more  than  10  minutes  without  severe  contamination  by  the  residual  gas  in 

the  vacuum  chamber.  Therefore,  all  the  ESD  experiments  discussed  in  this 

chapter  were  made  within  10  to  20  minutes  after  creating  a new  surface. 

Figure  19  shows  a desorption  rate  profile  of  m/e  = 16  (0+)  and  23 

(Na+)  signals  as  a function  of  electron  bombardment  time.  A 1 keV  scanning 

2 

electron  beam  with  a beam  current  density  of  1 mA/cm  was  used.  Note  that 
the  desorption  rate  for  m/e  = 23  (Na+)  rises  initially  for  the  first  20 
seconds.  This  is  followed  by  a monotomic  decrease  at  longer  electron  bom- 
bardment times.  For  m/e  = 16  (O"*") , the  initial  rise  was  not  observed. 
Instead,  the  rate  was  nearly  constant  for  20  seconds  and  then  decreased. 

The  absorption  current  did  not  vary  significantly  with  time  during  electron 
beam  bombardment,  suggesting  that  charging  of  the  specimen  did  not  occur. 

Similar  desorption  curves  were  obtained  for  m/e  = 23  (Na+)  for  energies 

“A  — ^ 

of  0.5,  1,  2 and  4 keV  beam  current  densities  varying  between  10  and  10 
2 

A/cm  . These  data  were  normalized  at  the  maximum  desorption  rate  for  each 
energy  and  are  plotted  versus  electron  dose  (the  product  of  beam  current 
density  and  electron  bombardment  time)  in  Figure  20(a),  (b),  (c)  and  (d) 


Primary  electron  beam  ener 
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Figure  19.  A desorption  rate  profile  of  m/e  = 16  and  23  signals  from  the  soda-silica  glass 

surface  as  a function  of  electron  beam  bombardment  time.  A 1 keV,  1 mA/cm^  electron 
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for  0.5,  1,  2 and  4 keV  electrons,  respectively.  Plotting  dose  on  the 

abscissa,  the  experimental  data  fall  on  a single  line  for  a fixed  energy 

-4  -3  2 

when  the  beam  current  was  varied  from  10  to  10  A/cm  . In  all  the 
results,  the  desorption  rate  curves  consist  of  an  initial  rise  followed 
by  a monotonic  decrease.  When  the  primary  beam  energy  was  increased,  the 
electron  dose  necessary  to  reach  the  maximum  rate  was  shifted  to  a higher 
value  and  the  rate  of  decay  after  the  maximum  was  decreased. 

The  maximum  desorption  rates  versus  beam  current  densities  between 

__  / __  o 2 

10  and  10  A/cm  (at  2 keV)  are  shown  in  Figure  21.  The  rate  becomes 

larger  as  the  current  density  increases,  but  the  relationship  is  not 

linear.  Similar  results  were  observed  at  other  primary  electron  energies. 

The  oxygen  desorption  signals  (m/e  = 16)  were  also  monitored  with 

-3  -2  2 

beam  current  densities  ranging  from  10  to  10  A/ cm  (at  2 keV) . The 
normalized  signals  versus  electron  dose  are  shown  in  Figure  22.  Again, 
most  of  the  data  form  a single  line  except  for  the  highest  beam  current 

_3  2 

density.  For  beam  current  densities  of  10  A/ cm  , the  desorption  rate 

-2  2 

does  not  vary  up  to  ~ 3 x 10  coulomb/cm  and  the  rate  subsequently 
becomes  lower,  similar  to  the  sodium  desorption  rate  at  2 keV.  The  oxy- 
gen KLL  Auger  signal  was  also  monitored  during  the  ESD  experiment.  The 
peak  height  decreased  with  increasing  irradiation  dose,  as  shown  in  Figure 
23.  The  oxygen  Auger  electron  escape  depth  is  at  least  several  times 
larger  than  the  ion  escape  depth  in  ESD  [66,88].  Therefore,  the  rate  of 
decrease  of  the  AES  signal  is  slower  than  for  the  ESD  signal. 

Although  the  data  show  that  the  electron  beam  desorbs  sodium  and 
oxygen  from  glass,  the  initial  rise  in  the  sodium  desorption  rate  cannot 
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Figure  20(a).  A desorption  rate  profile  of  m/e  = 23  (Na+)  signal  from 
the  soda-silica  glass  surface  as  a function  of  electron 
dose.  Figures  (a),  (b),  (c)  and  (d)  are  for  the  energies 
of  0.5,  1,  2 and  4 keV,  respectively,  with  varying  beam 
current  densities  shown  in  figures. 


Figure  20(b). 
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Figure  20(c). 
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Figure  20(d). 
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BEAM  CURRENT  DENSITY  (A/cm2) 


Figure  21.  A plot  of  maximum  desorption  rates  of  sodium  from  the 

soda- silica  glass  surface  as  a function  of  beam  current 
density  at  2 keV. 
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Figure  22.  A desorption  rate  profile  of  m/e  = 16  (0+)  signals  from  the  soda-silica 
glass  surface  as  a function  of  electron  dose  at  2 keV. 
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Figure  23.  A plot  of  oxygen  (KLL)  Auger  peak  height  from  the  soda-silica 
glass  surface  as  a function  of  electron  dose  at  2 keV. 
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be  explained  by  a conventional  ESD  mechanism,  where  the  rate  of  desorption 
is  expressed  by  an  exponential  decay  (equation  (2)).  Neither  the  variation 
of  rate  of  desorption  with  time,  nor  the  non-linearity  of  rate  with  beam 
current  density  can  be  explained.  This  is  because  both  electron  stimulated 
desorption  and  diffusion  processes  are  simultaneously  occurring  at  the 
surface.  A phenomenological  model  was  developed  and  solved  mathematically. 
The  model  is  compared  with  experimental  results  in  the  next  section. 

Phenomenological  Model  of  ESD  on  Glass  Surface 


A Model 

The  proposed  model  is  similar  to  one  used  by  Kirby  et  al.  [25, 

26]  to  explain  Auger  data  for  oxygen  or  carbon  monoxide  adsorbed  on  a 
silicon  surface,  where  the  adsorbed  molecular  species  are  dissociated 
under  the  electron  beam.  The  dissociated  species  may  then  desorb  or 
diffuse  into  the  silicon.  By  analogy,  the  following  process  is  considered. 
When  the  electrons  strike  the  glass  surface,  they  may  break  the  ionic  non- 
bridging oxygen-sodium  bonds  [98].  Ion  pairs  of  an  oxygen  ion  and  a sodium 
ion  are  created  at  the  surface.  The  oxygen  ion  traps  a hole  in  a non-bonding 
2p  orbital  [68r110]  and  the  sodium  is  no  longer  bound  to  the  oxygen  [68]. 

Such  defects  are  reported  in  studies  of  radiation  damage  in  soda-silica 
glass  [111].  Unless  there  is  some  mechanism  to  stabilize  the  defect,  the 
bond  will  reform  and  the  defect  will  disappear.  At  a surface,  however,  a 
stress  relaxation  of  broken  bonds  may  stabilize  the  defect  [112].  Thus, 
the  sodium  may  be  desorbed  into  the  vacuum  or  diffused  into  the  bulk. 

The  desorption  occurs  via  Coulombic  repulsion  with  0+,  which  is  produced 
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by  interatomic  Auger  de-excitation  or  intraatomic  Auger  de-excitation 
processes  [ 97  ,109]  (see  detailed  discussion  in  a later  section).  Sodium 
diffusion  proceeds  under  the  influence  of  the  electric  field. 


Mathematical  Formulation  of  Model 

The  processes  considered  above  are  described  in  the  following 
diagram. 

Na+  (desorption) 


I 

-Si-O-Na 

I 


-Si-0'Na+ 

l 


Na+  (diffusion) 


i 

where  -Si-O-Na  represents  the  state  in  which  the  sodium  and  non-bridging 
oxygen  is  ionically  bound,  and  -Si-0’  Na+  represents  the  state  in  which 
the  O-Na  bond  is  broken  by  electron  bombardment  and  a hole  is  trapped  at 

0. 

In  order  to  formulate  these  processes  mathematically,  the  following 
parameters  are  defined: 

k2  = a constant  associated  with  the  rate  of  conversion  of  sodium 
from  -Si-O-Na  to  -Si-0'  Na+. 

i i 

k^  = a constant  associated  with  the  rate  of  desorption. 

The  concentration  of  the  atoms  or  ions  is  time  dependent  as  defined  by: 
C^(t)  = surface  concentration  of  Na  bonded  to  non-bridging  oxygen 

i 

at  -Si-O-Na, 

i 


and 


61 


C0(t)  = surface  concentration  of  Na+  ion  after  conversion  at  -Si-0'Na+. 
1 i 

These  have  the  following  unit:  number  of  atoms  (including  ions)/cm^. 

If  we  assume  that  the  conversion  is  a first  order  reaction,  the 
rate  equations  for  the  system  are 


3“  Cx(t)  = -k^Ct)  (10) 

^ C2(t)  = -k3C2(t)  - ^ nb(t)  + k^t)  (11) 


where  the  number  of  atoms  which  have  diffused  into  the  bulk  is  given 
by  nb(t). 

In  evaluating  the  bulk  diffusion  term,  a one-dimensional  field-enhanced 
diffusion  equation  can  be  applied  as  discussed  in  Chapter  II. 


ft  VZ’C)  ‘ D ~2  Cb(2-C>  - F*  fz  V2^ 

3 2 


(12) 


where  F and  q are  electric  field  strength  and  mobility,  respectively. 

Here  z is  the  distance  perpendicular  to  the  surface.  A quantity,  Cb(z,t), 
is  the  bulk  concentration  with  a dimension  of  number  of  atoms/cm^.  A 
parameter,  p,  with  a dimension  of  cm  ^ is  introduced,  where  p is  a distance 
at  which  sodium  is  considered  to  be  bulk  sodium,  as  opposed  to  surface 
sodium  available  for  ESD  or  Auger  measurement.  The  boundary  condition  for 
equation  (12)  is  taken  as 


Cb(0,t)  = p C2(t)  . 

The  intital  conditions  are 


(13) 


Cb(0)  = C 


01 


and  C2(0)  = 0 


(14) 
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where  is  the  number  of  sodium  atoms  bonded  to  non-bridging  oxygen 
before  electron  bombardment. 

In  Appendix  B,  we  obtain  the  solution  of  equation  (12)  as 


C (z  t)  = pz  • exp(^  z - t) 

b 2/rrD  2D  4D 


(15) 


X 


f C„  (t')  • exp((Fy)2t' /4D)  • (t-t')  3^2  • exp(-z2/4D(t-t'))dt' . 

0 

The  number  of  atoms  in  the  bulk  (semi-inf inite  solid)  is  then  given  by 

03 

n,  (t)  = f C (z,t)dz  (16) 

b 0 b 


• exp(-^M^t).  / c2(t’)  • exp(i|g^t’).  (t-t’)‘3/2 


2/rrD 
• Q(t , t ' )dt ' 


(17) 


where 


t 1/2 

Q(t , t ' ) = 2D(t-t’)  {l+(-^ili Fy  • exp(^|^(t-t')) 


/D 


Erf c (— ^—  (t-t’)1/2)}  . 

2/D 


Differentiating  with  respect  to  t,  we  obtain 


dnb(t) 

dt 


p/D  * 


,.x-3/2 


(Fy)‘ 


(18) 


- - * / C9  (t ' ) • (t-t  * ) • exp(--^-(t-t’))dt\  (19) 


2/7  0 


Substituting  (19)  and  the  solution  of  (10)  into  (11) , we  obtain  the 
coupled  integro-dif ferential  equation  for  the  model,  as  solved  below. 
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Solution 

The  exponential  term  in  equation  (19)  may  be  expanded  by  a series 
approximation  and  we  integrate  by  parts  using  the  conditions  of  (13) 
and  (14).  We  obtain 


dnb(t) 

dt 


p/D 

2/ TT 


E 

n=0 


t dC2(t’) 


dt' 


(-Dn  / (Fy) 2.  n (t-t')n  1/2 
n!  ^ 4D  ; n-1/2 


dt'  . 


(20) 


Inserting  (20)  into  (11),  the  Laplace  transform  of  equation  (11)  is 
obtained  as 


s C2(s)  = -k3C2(s) 


L{ 


dnb(t) 

dt 


} + 


k2C01 
s + k„ 


(21) 


where  C2(s)  is  the  Laplace  tranformation  of  the  concentration  C2(t),  i.e.. 


C2(s)  = £{C2(t)}. 


(22) 


Using  the  convolution  theorem  in  the  Laplace  transformation,  equation 
(21)  is  reduced  to  (Appendix  C) , 


s C2(s)  - -k^C^s)  + 


" ,(-l)n  ,(Fp)2.n  1 

1 ■ J7T72 


p/D 

2/u"  n=0  nT 


. S (s)}+  ^£01 

n-1/2  L2(s;j+s  + k9 
s 2 

where  T(n+l/2)  is  the  Gamma  function. 

Solving  equation  (23)  for  C2(s),  we  obtain 

)/D* 


(s)  = k2C01(s+k2)’1*[s'Ht^E“  * E {'T^<1IS2-)n*^T77 

(24) 


2 2 01v  2/  "3  2i/rjj-  n^QL  n!  v 4D  7 n-1/2  _ n-1/2 


(23) 


If  we  replace  the  infinite  summation  in  equation  (24)  by  a partial 
summation,  which  is  taken  large  enough  to  minimize  the  error,  we  obtain 
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c2(s)  * k2C01(s+k2)-1.  [3*3-  . £ {lliil1 . 

2/  tt  n=0 


r(n+l/2)  j ]-l 


n-1/2  n-1/2 

s 


We  can  then  write 


K 


C2(s)  = 


k2C01s 
/s"  M(s) 


where 


M(s)  - (S*2).  (SK+1/2  + ^ + £ „ . 3K-")  . 

1 3 A 11 

n=0 

The  expression  for  M(s)  can  be  written 
K+3 

M(s)  - II  (/s  + g.) 
j=l  3 

where  the  g are  constants  to  be  determined.  Using  equation  (27) 
rearranging,  we  can  write  equation  (26)  as 
K+3  B . 

C2(s)  = Z -J 

j=l  /^(/ibg^) 

where  the  are  constants  to  be  determined.  The  inverse  Laplace 
form  of  each  term  is  given  by  [113] 


L 1 “ = F( g/t) 

/s(/s+g) 


where 


F(g/t)  = exp(g2t)  • erf c (g/t) . 


(25) 


(26) 


(27) 


(28) 

and 


(29) 

trans- 


(30) 


(31) 
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It  is  difficult  to  use  this  solution  (equation  (29))  since  we  need  to 
solve  equation  (27)  for  M(s)  and  obtain  the  roots,  g^  , then  solve  for  B_. 
from  equation  (29).  The  choice  of  K is  largely  dependent  on  the  magni- 
tude of  (Fy)2/4D  in  equation  (25).  From  our  knowledge  of  previous  works, 
(Fy)  /4D  is  in  the  order  of  10  near  room  temperature.  Therefore,  K=2 
was  chosen  as  a first  order  approximation.  Equation  (26)  is  then  reduced 
to 

2 

C?(s)  = /s"  k2CQ2*  (s+k2)  {s /s  + p i/5~ * s + k^/s”  + p/D*(^gj^  )}  . (32) 

The  roots,  g.  and  B.  are  obtained  in  the  same  manner  as  discussed  above. 

2 2 

The  time  dependence  of  the  concentration  C2(t)  can  be  written  in  a 
reduced  manner  as 


c2(t) 


5 

2 B f(  g./t). 

j=l  3 3 


(33) 


The  mass  spectrometer  signal  is  directly  proportional  to  the  desorption 
rate,  and  may  be  expressed  as 


IESD^  " Bk3*C2(t) 

where  6 is  a proportionality  constant. 


(34) 


Model  Calculation 

A quantitative  comparison  of  this  model  with  the  experimental  data 
is  difficult,  since  very  little  is  known  about  the  rate  coefficients,  k2 
and  k^.  We  have  previously  used  literature  data  to  estimate  the  diffusion 
coefficient  of  sodium  in  the  glass  and  determined  the  magnitude  of  the 
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electric  field  induced  by  the  electron  beam  (Chapter  II) . Using  these 
values,  the  rate  coefficients,  k ^ and  k^,  can  be  evaluated  by  a best  fit 
to  the  experimental  data. 

In  general,  the  rate  coefficient  is  dependent  upon  the  beam  current 
density  and  energy.  For  a fixed  energy,  it  is  expressed  [70]  as 

k = jpo/e  (35) 

where  j and  e are  the  beam  current  density  and  the  electronic  charge, 

respectively,  and  a is  the  conversion  or  desorption  cross  sections.  The 

cross  section  is  a function  of  energy,  thus  the  conversion  and  desorption 

rates,  and  k^,  are  respectively  expressed  as  the  conversion  cross 

section,  a , and  the  desorption  cross  section,  a„OT..  As  previously 
conv.  ESD 

-15  2 5 

estimated,  the  value  of  D and  F were  taken  as  1 x 10  cm  /sec  and  1 x 10 

V/cm  and  the  mobility  is  calculated  by  the  Einstein  relation,  as  discussed 

in  Chapter  II.  The  parameter,  p,  was  equated  to  the  reciprocal  of  the  ESD 

escape  depth  and  was  estimated  to  be  - 4 x 10 ^ cm  ^ since  ESD  originates 

only  from  the  surface  monolayer  [88].  The  constants,  k 2 and  k^,  were 

adjusted  for  a best  fit  to  the  experimental  data  and  the  theoretical 

desorption  rate  versus  the  electron  bombardment  time  was  calculated.  This 

is  shown  in  Figure  24,  where  the  calculated  curves  were  normalized  to  their 

maximum  rates.  The  calculated  curves  are  compared  to  experimental  data 

2 

corresponding  to  a beam  current  density  of  1 mA/cm  at  various  primary 
electron  energies  (Figures  18(a),  (b) , (c)  and  (d)).  The  error  is  attri- 
buted to  statistical  variations,  but  they  are  systematic  deviations  when 
considering  a long  period  of  time.  The  time  required  to  reach  the  maximum 


67 


OJ 

o 

cO 

y-i  • 
M O 
3 CM 
CO 

0) 

CO  M 
CO  3 
CO  00 

t-H  *H 

00  Fn 

co  e 
a o 

•H  h 


co  -a 
I 3 
3 a 
T3  3 
O T3 
CO  O 


<1) 

J3 


t-i 

a 

3 

M 

CO 


e 

o 
)-i  co 

M— I "O 

§ 


cd 
u 
■a  c 
o a 
cn  s 

•H 
■W  J-i 
O Q) 
CL 

cn  x 
a)  <u 


cd  a; 
M s: 
■u 

c 

o o 


C-  -U 

O M-l 

cn 

O)  4J 

'td  cn 
aj 

T3  -Q 

aj 

U QJ 

u 
cd 
3 

a 4J 

rH  Cd 

«d  si 

CJ 


cd 


CN 

0) 

P 

60 


68 


desorption  rate  and  the  decay  rate  are  well  matched  by  the  calculated 

data.  The  conversion  cross  section  and  the  desorption  cross  section  for 

each  energy  are  evaluated  from  equation  (35)  and  are  summarized  in  Table 

“18  2 

1.  The  cross  sections,  a and  o_„,  are  on  the  order  of  10  cm 

conv.  ESD 

-19  2 

and  10  cm  , respectively,  and  decrease  as  the  energy  increases. 

Mechanism  For  Sodium  Desorption  From  Glass  Surface 
In  the  previous  sections,  it  has  been  shown  that  sodium  and  oxygen 
are  desorbed  from  glass  surfaces  and  the  kinetics  have  been  modeled  in 
a phenomenological  manner.  However,  the  mechanism  of  sodium  desorption 
is  not  very  obvious.  With  the  conventional  Menzel-Gomer-Redhead  ESD 
model  [ 95  , 96  ] sodium  desorption  from  the  soda-silica  glass  cannot  be 
fully  explained.  Instead,  a mechanism  proposed  by  Knotek  and  Feibelman 
[ 97  ,109],  an  interatomic  Auger  de-excitation  model,  can  be  applied  to  the 
present  case. 

Knotek  and  Feibelman  [ 97  ,109]  have  formulated  a mechanism  for  ESD 

of  positive  ions  from  maximal  valence  transition  metal  oxides,  i.e., 

oxides  in  which  the  metal  cation  is  ionized  to  the  rare  gas  configuration 
4+ 

(e.g.,  Ti  in  Ti02).  In  the  metal  oxides  with  a high  negative  charge  on 
the  oxygen,  the  production  of  positive  oxygen  is  only  possible  by  multiple 
ionization.  An  interatomic  Auger  process  can  cause  multiple  ionization, 
starting  with  the  creation  of  an  inner  level  electron  hole  on  the  cation 
which  is  filled  by  an  Auger  decay  involving  the  oxygen  valence  electrons. 
This  leads  to  a loss  of  two  electrons  on  the  oxygen  site.  An  Auger  cas- 
cade can  lead  to  an  even  higher  positive  charge.  A schematic  illustration 


69 


Table  1.  Summary  of  Calculated  Results  for  the 
Conversion  and  the  Electron  Stimulated  Desorption  Cross  Sections 


Eo 

500  eV 

1000  eV 

2000 

eV 

4000  eV 

C\| 

0.09 

0.05 

0.04 

0.025 

k3 

0.0065 

0.005 

0.003 

0.002 

a 

conv. 

1.4  x 10-17 

8.0  x 10"18 

6.4  x 

IQ'18 

4.0  x 10~18 

°ESD 

1.0  x 10  18  8.0  x 10~19 

Unit  of  cross 

4.8  x 
section  ( 

io-19 

cm2) 

-19 

3.2  x 10 
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of  an  interatomic  Auger  decay  in  Ti02,  for  example,  is  shown  in  Figure 
25.  The  highest  occupied  level  of  the  Ti^+  is  Ti-3p  level  of  ~ 34  eV 
since  there  are  no  electrons  in  the  valence  band.  If  an  electron  ionizes 
this  level,  the  dominant  decay  mode  will  be  an  interatomic  Auger  decay, 
since  the  valence  electron  density  on  the  Ti  ion  is  too  low  for  intra- 
atomic  Auger  decay  to  be  a likely  process.  Therefore,  one  of  the  0-2p 
electrons  decays  into  the  Ti-3p  hole  and  another  electron  is  ejected 
from  the  0-2p  level  in  the  Auger  process.  The  0“  ion  is  thus  converted 
to  an  0 . If  one  of  the  remaining  valence  electrons  is  shared  between 
the  0 and  Ti  (recall  that  in  equilibrium  an  0 is  not  02_  but  0^1+X^~) 
then  there  is  an  appreciable  probability  that  0 will  desorb  as  0+. 

This  picture  is  now  applied  to  our  case.  When  the  freshly  created 
soda-silica  glass  surface  is  bombarded  by  electrons,  they  will  break  the 
O-Na  bond  and  the  unsatisfied  non-bridging  oxygen  bond  will  trap  a hole 
in  a non-bonding  2p  orbital  of  oxygen,  resulting  in  a Si4+01_Na1+  configur- 
ation. If  the  sodium  desorbs  as  Na4+,  it  only  occurs  via  Coulombic  repul- 
sion with  0 ions.  Therefore,  oxygen  has  to  be  converted  into  positive 
oxygen.  This  can  occur  either  by  intraatomic  Auger  de-excitation  at 
oxygen  or  interatomic  Auger  de-excitation  as  applied  to  Si4+01-.  An 
energy  diagram  for  Si4+01_Na1+  is  given  in  Figure  26. 

Intraatomic  Auger  de-excitation.  The  highest  occupied  level  of  Cf 
is  0-2p  level,  which  forms  a valence  band.  From  photoemission  measurements 
of  the  valence  level  of  amorphous  SiC^  [114]  the  valence  band  consists  of 
non-bonding  and  bonding  bands.  These  are  considerably  broad  with  band 
widths  of  approximately  3.2  eV  and  5.7  eV,  respectively.  When  the  incident 
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Figure  25.  Interatomic  Auger  de-excitation  model  for  stimulated 

desorption  of  ions.  Formulation  of  a core  hole  in  the 
3p  shell  of  maximal  valency  Ti02  is  followed  by  inter- 
atomic Auger  decay. 
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Figure  26.  Interatomic  Auger  de-excitation  model  for  stimulated 
desorption  of  sodium  ions.  After  a hole  trapping  in 
non-bridging  oxygen  2p  orbital  of  a single  oxygen, 
formulation  of  a core  hole  in  the  2p  shell  of  maximal 
valency  of  Si  is  followed  by  interatomic  Auger  decay. 
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electron  ionizes  an  0-2s  level  (L^) , an  electron  from  either  the  bonding 
or  the  non-bonding  band  occupies  the  0-2s  shell  hole  and  the  excess  energy 
is  released  in  the  de-excitation  by  emission  of  another  electron  from  the 
valence  band  (O-LW  Auger  electron) . The  oxygen  0 is  thus  converted  to 


Interatomic  Auger  de-excitation.  The  core  level,  Si-2p(I,2  of 
- 100  eV  is  the  highest  occupied  level  in  Si  since  no  electrons  are  in 
the  valence  band.  If  an  incoming  electron  ionizes  this  level,  one  of  the 
electrons  in  the  0 valence  band  (e.g.,  a 2p  non-bonding  electron)  drops 
to  occupy  the  Si-2p  level  by  interatomic  transition.  The  energy  released 
by  this  process  is  approximately  -100  - -3  = -97  eV,  which  is  high  enough 
to  eject  another  electron  (0-L?^_V^V^  Auger  electron).  This  multiple 

Z , J 

ionization  process  leaves  oxygen  as  0+. 

The  positive  oxygen  created  by  either  intraatomic  or  interatomic 

Auger  de-excitation  processes  results  in  Coulombic  repulsion  with  Na+. 

Thus,  the  sodium  desorbs  as  Na+.  In  addition  to  sodium  desorption,  oxygen 

4-  4+ 

has  to  be  desorbed  as  Cr  via  Coulombic  repulsion  with  Si  . Desorption 
of  0+  was  experimentally  observed,  as  discussed  in  the  previous  section. 

This  hypothesis  is  supported  by  the  data  on  Na+  desorption  rate 
versus  the  primary  electron  beam  energy,  as  shown  in  Figure  27.  The 
sodium  desorption  rate  above  - 100  eV  increases  rapidly.  This  indicates 
that  the  sodium  desorbs  mainly  via  the  interatomic  Auger  de-excitation 
process  since  the  ionization  threshold  of  the  Si-2p  level  (-  100  eV) 
coincides  with  the  primary  electron  beam  energy  above  which  the  sodium 
desorption  rate  increases  drastically.  Below  100  eV,  a small  sodium 
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Figure  27. 


Sodium  desorption  yield  from  a soda-silica 
glass  surface  versus  electron  energy. 


75 


rate  was  observed  but  the  signal  intensity  was  weak.  Below  60  eV  it  was 
too  small  to  be  measured.  Hence,  the  threshold  energy  for  intraatomic 
Auger  transition,  which  should  be  equal  to  the  0-2s  ionization 

energy,  20.2  eV,  could  not  be  found  experimentally.  This  data  indicates 
that  sodium  desorption  occurs  via  an  interatomic  Auger  de-excitation 
process. 


Discussion 

As  shown  in  the  previous  calculation,  we  have  seen  a reasonable 
agreement  between  the  experimental  data  and  the  model  calculations.  A 
conventional  ESD  rate  equation  was  not  adequate  to  explain  the  sodium 
desorption  kinetics,  especially  the  initial  rise  observed  in  the  sodium 
desorption  rate.  An  obvious  question  is  whether  the  initial  rise  in 
kinetics  is  due  to  surface  charging  or  contamination.  For  example,  the 
absorption  current  measured  between  ground  and  the  specimen  during  the 
ESD  experiment  reached  a steady  value  within  a few  seconds,  while  the 
desorption  rate  maximum  was  observed  after  10  to  50  seconds,  depending 
on  the  beam  current  density.  Further  bombardment  did  not  vary  the  absorp- 
tion current  significantly.  In  addition,  the  initial  rise  was  much 
smaller  for  oxygen  desorption  rate.  In  fact,  the  rate  was  constant  over 
the  period  which  the  sodium  desorption  rate  was  increasing.  If  the  ini- 
tial rise  is  due  to  the  time-dependent  surface  charging,  the  absorption 
current  should  have  varied  [115]  and  the  oxygen  desorption  rate  should 
have  increased,  as  did  the  sodium  desorption  rate.  Furthermore,  the  Auger 
peak  energy  did  not  vary  during  electron  bombardment.  This  further  substan- 
tiates that  the  surface  charge  was  not  time-dependent,  but  only  stable 
charging  occurred. 
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It  may  also  be  argued  that  the  initial  rise  results  from  surface 

contamination.  As  shown  in  Figure  18,  the  surface  was  covered  with 

-9 

F and  CO  after  90  minutes  in  a vacuum  of  1 x 10  Torr.  The  flux  of 
residual  gas  molecules  which  strike  the  surface  per  square  centimeter 

-9 

per  second  is  proportional  to  the  gas  partial  pressure  [70].  At  1 x 10 

11  2 

Torr  the  flux  of  gas  molecules  is  about  10  molecules/cm  sec.  There- 

14 

fore,  after  90  minutes,  approximately  5.4  x 10  molecules  of  residual 

gas  have  collided  with  one  square  centimeter  of  the  fracture  surface. 

15  2 

Since  a monolayer  is  about  10  molecules /cm  , a 90-minute  exposure  to 
-9 

1 x 10  Torr  will  result  in  less  than  a monolayer,  even  if  the  sticking 

coefficient  is  unity  for  all  molecules.  The  ESD  experiment  was  carried 

out  within  10  to  20  minutes  after  creation  of  the  surface.  Therefore, 

the  contamination  on  the  surface  is  less  than  one-tenth  of  a monolayer 

for  a sticking  coefficent  of  unity.  The  ESD  cross  section  for  adsorbate 
“16  2 

is  about  10  cm  [88].  With  such  a cross  section,  63%  of  the  adsorbates 
would  be  desorbed  in  one  second  of  electron  bombardment.  Thus,  the  conta- 
mination is  not  a problem  in  this  case,  nor  can  it  be  responsible  for  the 
initial  rise  in  the  sodium  desorption  rate. 

Mechanisms  of  sodium  ESD  from  glass  surfaces  were  studied  and  two 
processes  were  considered  in  order  to  create  0'''+  from  0^  . These  were 
intraatomic  Auger  de-excitation  at  oxygen  and  interatomic  Auger  de-exci- 
tation between  Si  and  0.  From  the  threshold  measurement,  it  was  found 
that  the  interatomic  Auger  de-excitation  was  the  most  likely  process. 
However,  there  remains  a question  as  to  why  the  intraatomic  Auger  decay 
of  the  oxygen  core  hole,  which  also  leads  to  the  removal  of  two  valence 
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electrons  and  thus  to  the  same  final  state  as  the  Si-2p  hole  decays, 
should  not  lead  to  a stronger  threshold  than  the  interatomic  Auger  de- 
excitation process.  This  can  be  explained  by  the  law  of  conservation 
of  energy.  In  the  intraatomic  Auger  de-excitation  process,  the  maximum 
energy  released  by  the  2s  hole  is  20.2  - 3.0  (energy  of  the  valence  band 
maximum)  = -17.2  eV  (Figure  26).  In  order  to  provide  the  maximum  energy 
to  the  desorbing  atom  (oxygen) , we  suppose  that  the  second  electron  in 
the  Auger  process  is  kicked  into  the  conduction  band,  which  has  the  mini- 
mum energy.  In  the  absence  of  "hole-hole  interaction,"  this  energy  would 
be  the  band  gap  energy,  i.e.,  ~3  eV  above  the  Fermi  energy  level  [114]. 
However,  it  is  known  that  a second  electron  from  the  oxygen  sitting  in 
the  Madelung  potential  will  require  more'  energy  than  that  to  remove  the 
first  valence  electron.  This  additional  energy  is  just  the  "hole-hole 
repulsion"  energy.  This  may  be  evaluated  by  a comparison  of  KW  Auger 
data  with  x-ray  photoemission  results,  to  be  on  the  order  of  6 eV  [116]. 
Thus,  the  minimum  energy  required  to  remove  the  second  valence  electron 
from  an  0 atom  is  then  6+3+6  =15  eV.  Therefore,  the  kinetic  energy 
of  the  0+  is  only  17.2  - 15  = 2.2  eV.  By  contrast,  much  more  kinetic 
energy  is  available  in  the  interatomic  Auger  process,  since  the  Si-2p 
ionization  energy  is  about  100  eV,  which  is  - 80  eV  higher  in  the  kinetic 
energy  of  desorbing  0+  ions.  Consequently,  the  desorption  via  interatomic 
Auger  de-excitation  is  the  most  likely  process. 

Because  of  the  large  charge  transfer  involved  in  the  ESD  process,  it 
may  be  anticipated  that  ESD  data  provides  further  information  bearing  on 
the  states  and  valencies  of  bonding  in  the  outermost  layer  of  a solid. 
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This  may  be  extremely  useful  to  study  the  early  stage  glass  surface  reaction 
with  environment  (for  example,  weathering  and  ion  exchange  processes). 

This  unique  capability  of  ESD  is  a consequence  of  the  facts  that: 

(1)  a core  electron  must  be  excited  to  initiate  desorption,  thus,  giving 
rise  to  a bonding-state-dependent  threshold  energy,  and  (2)  a nearest 
neighbor  atom  is  involved  in  an  interatomic  Auger  de-excitation,  which 
results  in  valence  sensitivity.  Thus,  ESD  can  be  a probe  of  the  valence 
of  bonding  at  a surface  and  it  is  promised  that  careful  threshold  measure- 
ment will  provide  important  information  directly. 

The  ESD  cross  section  can  be  estimated  based  on  the  fundamental 
mechanism.  The  interatomic  Auger  process  is  one  in  which  a primary  elec- 
tron creates  a core  hole  of  a parent  atom  followed  by  an  Auger  neutrali- 
zation from  a nearest  neighbor  atom.  Therefore,  the  sodium  ESD  cross 
section  is  equal  to  the  product  of  the  inner  shell  ionization  cross 
section  of  Si-2p  and  the  probability  of  sodium  via  Coulombic  repulsion 
is  unity.  McGuire  [117]  has  calculated  the  Si-2p  ionization  cross  sec- 
tion based  upon  the  Bethe  theory  [ 99  ] . For  primary  energies  between 

■”18  ”*18 

0.5  and  4.0  keV,  ionization  cross  sections  of  - 4 x 10  to  ~ 2 x 10 
2 

cm  were  reported.  A similar  value  of  cross  section  was  obtained  experi- 
mentally by  Meyer  and  Vrakking  [118,119].  Since  the  interatomic  transi- 

-14 

tion  rate  is  very  fast  (-  10  seconds  [120]),  the  probability  of  the 
Auger  neutralization  is  only  limited  by  the  density  of  state  of  the 
oxygen  valence  band.  Because  of  the  localized  electron  states  in  the 
oxygen  valence  band,  i.e.,  non-bonding  and  bonding  electrons,  the  proba- 
bility of  the  Auger  neutralization  is  low  [121].  Furthermore,  the  survival 
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probability  for  0+  via  Auger  neutralization  is  - 0.01  to  ~ 0.1  [120]. 

-19  2 

Thus,  the  estimated  sodium  ESD  cross  section  is  ~ 10  cm  or  less. 

This  is  still  smaller  than  the*  value  determined -by '-the  experi- 
-19  2 

ment  (~  10  cm  ) ; however,  the  calculation  provides  a lower  limit  on 
the  possible  ESD  cross  section. 

In  the  ESD  threshold  measurement , the  energy  was  only  scanned  up 

to  200  eV,  therefore,  only  the  oxygen  L 2 transition  was  possible. 

If  the  energy  is  increased,  the  oxygen  KLL  transition  would  occur.  It 

is  speculated  that  the  desorption  rate  would  increase  more  when 

the  oxygen  K edge  was  crossed  at  530  eV.  Because  of  this,  the  probability 

for  oxygen  multiple  ionization  would  increase  with  higher  energy,  resulting 

in  a higher  sodium  desorption.  Therefore,  the  observed  sodium  ESD  cross 
-19  2 

section  of  10  cm  is  very  reasonable. 

In  electron  stimulated  desorption  of  sodium,  the  desorption  rate 

increased  when  the  beam  current  density  was  increased  (Figure  21) , but 

not  linearly.  Based  upon  the  present  model,  the  dependence  of  the  maximum 

desorption  rate  upon  current  density  is  evaluated  from  equation  (19)  and 

is  shown  in  Figure  28.  For  a small  beam  current  density,  the  desorption 

rate  was  negligible,  and  thus,  ESD  is  not  important.  As  the  current 

-3  2 

density  is  increase  to  more  than  10  A/ cm  , the  maximum  desorption  rate 

increases.  The  desorption  rate  does  not  continue  to  increase,  but  tends 

to  saturate  at  a very  high  beam  current  density.  This  means  that  ESD 

-3  2 

becomes  important  at  beam  current  density  > 10  A/cm  , but  at  beam 

-2  2 

current  densities  ~ 10  A/cm  diffusion  limits  further  increase  in  the 
ESD  rate.  In  a typical  AES  experiment,  a beam  current  density  range  of 
10  to  10  A/cin  is  used.  This  is  the  range  where  ESD  becomes  important. 
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Electron  Beam  Current  Density  (A/cm2) 


Figure  28.  Calculated  maximum  desorption  rate  of  sodium  from  the 

soda-silica  glass  surface  as  a function  of  beam  current 
density. 
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While  the  proposed  model  correctly  predicts  the  rise  and  the  fall 
in  the  sodium  desorption  rate,  it  predicts  that  the  rate  should  be  zero 
at  t = 0 (Figure  24) . The  initial  conditions  we  have  used  to  obtain 
the  solutions  are  applied  to  ideal  surface  conditions,  i.e.,  no  defects 
at  the  surface  before  electron  bombardment.  However,  this  may  not  be 
true  in  a real  case.  Fracturing  is  essentially  a bond  break  and  defects 
are  formed.  Therefore,  it  is  probable  that  non-bonded  sodium  exists 
prior  to  the  electron  bombardment.  In  this  case,  the  initial  condition 
(14)  should  be  changed  to 

C2(0)  = C02  (36) 

where  is  the  concentration  of  non-bonded  sodium  ions  before  electron 

bombardment.  However,  solving  the  equation  with  this  initial  condition 
is  not  easy  because  a discontinuity  how  exists  at  z = 0,  making  it  diffi- 
cult to  find  the  Laplace  transform  of  C^(z,t).  This  is  one  difficulty. 

In  addition,  the  concentration  of  non-bonded  sodium  ions,  is  diffi- 

cult to  quantify  and  varies  from  one  experiment  to  the  next.  Because  of 
this,  a calculation  using  the  appropriate  initial  conditions  was  not 
attempted. 


Summary 

1.  Electron  stimulated  desorption  of  sodium  from  in  situ  fractured  soda- 
silica  glass  surfaces  has  been  studied.  The  desorption  rate  was 
dependent  upon  the  electron  beam  energy,  current  density  and  bombard- 


ment time. 
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2.  The  desorption  rate  of  sodium  first  increased,  followed  by  a mono- 
tonic decrease.  This  initial  rise  in  the  rate  cannot  be  explained 
by  a conventional  ESD  rate  equation.  It  was  found  that  the  initial 
rise  in  the  desorption  rate  was  an  inherent  characteristic  of  ESD 
from  glass. 

3.  A phenomenological  model  was  developed  considering  processes  involving 
production  of  ion  pair,  followed  by  a hole  trap  at  the  oxygen  site 
and  electron  stimulated  desorption  and  field-enhanced  diffusion.  The 
model  predicts  correctly  the  kinetics  of  sodium  desorption. 

4.  The  electron  stimulated  desorption  cross  section  was  evaluated  to  be 

-19  2 

about  10  cm  by  comparing  the  experimental  data  with  the  model. 

It  decreased  when  the  energy  was  increased. 

5.  A mechanism  for  sodium  desorption  from  soda-silica  glass  was  developed 

by  considering  intraatomic  and  interatomic  Auger  de-excitation  processes. 
For  the  primary  electron  energies  above  100  eV,  the  interatomic  Auger 
de-excitation  process  is  the  most  likely.  ESD  cross  section  was  esti- 
mated based  upon  the  fundamental  principle  and  it  was  shown  that  the 
value  of  experimentally  determined  ESD  cross  sections  was  very  reason- 
able. 

6.  Because  of  the  large  charge  transfer  involved  in  the  ESD  process,  ESD 
can  be  a valence-sensitive  probe  at  a surface.  Careful  threshold 
measurement  may  provide  important  information  of  a surface  directly. 

7.  ESD  is  important  in  a conventional  AES  analysis  when  beam  current 

-3  2 

densities  > 10  A/cm  are  used. 


CHAPTER  IV 

THEORETICAL  INTERPRETATION 
OF 

ALKALI  AUGER  SIGNAL  DECAY  PROCESS 


We  have  shown  that  the  decay  of  an  alkali  Auger  signal  at  the  glass 
surface  during  analysis  with  Auger  electron  spectroscopy  was  a result  of 
its  diffusion  into  the  bulk  and  its  desorption  from  the  surface  under  the 
influence  of  the  electron  beam.  The  magnitudes  of  these  two  effects  were 
dependent  upon  electron  beam  energy,  beam  current  density,  substrate  tem- 
perature and  amount  of  contamination  on  the  surface.  These  parameters 
have  to  be  carefully  controlled  for  accurate  and  precise  Auger  analysis. 

The  observation  of  an  "incubation  time"  during  AES  of  the  glass  sur- 
face [39]  is  interesting.  During  the  incubation  time,  the  alkali  concen- 
tration can  be  determined  if  the  AES  measurement  is  accomplished  within 
this  period.  The  length  of  the  incubation  time  varies  strongly  with  the 
electron  beam  conditions,  the  glass  composition-  and  the  temperature. 

These  features  have  been  found  experimentally.  However-,  neither  a compre- 
hensive explanation  of  the  origin  of  the  incubation  time  nor  a theoretical 
prediction  has  been  made. 

In  this  chapter,  we  will  develop  a theory  describing  the  alkali  Auger 
signal  decay  process.  The  theory  will  be  applied  to  low  temperature  measure- 
ments and  the  origin  of  the  incubation  time  will  become  apparent. 
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Mathematical  Formulation 
of 

The  Alkali  Auger  Signal  Intensity 
In  Auger  electron  spectroscopy,  the  Auger  electrons  from  atoms  on 
the  surface  and  within  the  outermost  atomic  layers  of  the  surface  can 
escape  without  collisions  and  contribute  to  the  Auger  signal.  We  assume 
that  the  Auger  intensity  is  proportional  to  the  concentration  of  the 
atoms  on  the  surface  and  atoms  near  the  surface  which  have  diffused  into 
the  bulk.  Hence,  from  the  model  discussed  in  Chapter  III,  the  intensity 
of  the  alkali  Auger  signal  is 

1 (t)  Auger  = “^4(1)  + C2(t)  + nA(t)]  (37) 

A 

where  a is  a proportionality  constant  and  n (t)  is  the  number  of  bulk 
atoms  whose  Auger  electrons  escape  from  the  surface  without  collision. 

This  latter  quantity  is  expressed  as 

A °° 

n (t)  = / P(z)  • Cb(z,t)  ♦ dz  (38) 

o 

where  P(z)  is  the  probability  that  the  Auger  electron  will  escape.  We 
take  [122] 

P(z)  = exp(-z/XA)  (39) 

O 

where  X is  the  Auger  escape  depth  (on  the  order  of  10  to  50  A [66]). 
Characterization  by  this  simple  form  is  an  idealization  of  the  probability 
which  depends  on  the  details  of  the  scattering  process,  but  this  approxi- 
mation has  been  shown  to  work  well  [122].  The  concentration  of  atoms  in 

A 

the  bulk  is  given  by  equation  (15),  derived  in  Chapter  IV,  thus,  n (t)  is 


expressed  as 
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nA(t)  = / / exp(-z/A  ) • • exp(^  z - • t)  • C2(t') 

0 0 A 2^D  ZU 

2 2 

X exp(^|^  t')  *(t-t')'3/2-exp(-  |5(t-t,)_1)dzdt'.  (40) 

In  order  to  evaluate  the  integral,  we  put 

m = 1/X  - Fy/2D  (41) 

and 

q = D(t-t')  (42) 

and  we  use  the  fact  that  [113] 

°o  2 2 

f exp  (-mz)*  z *exp  (-z /4q)dz  = 2q{-/irq  "Hi'e^111  • erfc(/q  m)+l).  (43) 

0 

One  can  perform  the  integration  over  z in  equation  (40),  resulting  in 

nA(t)  = -2-  -h(t)  (44) 

/ttD 

where 

h(t)  = D f {1- /T-yF(y)  } • C,(t')-exp(-  (t-t ' ) ) • (t-t ' )_1/2dt ' . 

0 (45) 

In  equation  (45) 

y = m ]/q~  (46) 


and 


V2 

F(y)  = ey  • erf(y). 


(47) 
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We  restrict  our  calculation  of  the  Auger  intensities  to  a small  amount 
of  time.  When  considering  only  at  short  times,  we  can  reduce  equation 
(44)  to 


nA(t)  - ^ j C (t')-(t-t'  r1/2.exp(-  ^-(t-t’))dt\ 


(48) 


^ 0 


Now  consider  a case  where  a specimen  at  liquid  nitrogen  temperature 
(110 °K)  was  analyzed.  The  estimated  diffusion  coefficient  of  sodium  is 

o . , ,,-28  2 i-j 

about  10  cm  /sec  and  the  mobility  at  this  temperature  is  ~ 10  cm  /V 
sec.  Assuming  that  the  electric  field  is  not  influenced  by  the  tempera- 

c 2 

ture,  we  take  F = 10  V/cm.  With  these  values,  (Fy)  /4D  is  calculated 
to  be  - 10-17  and  the  exponential  term  which  appeared  in  equation  (48) 
does  not  vary  at  short  times.  Therefore,  we  can  put 


exp (- (Fy ) (t-t ' ) / 4D)  = 1. 


(49) 


Equation  (48)  can  be  reduced  further  to 


n.  = p 

low  temp. 


b -I/? 

= p/D/A  f c0(t')*(t-t’)  x/  dt'. 


(50) 


0 


Recalling  the  convolution  theorem  described  in  Appendix  C,  we  can  rewrite 
equation  (50)  as 


nA  „ * p/D//7-L  ^{D(C?(t))-L(t  1/2)}. 

low  temp.  l 


(51) 


The  Laplace  transform  of  C„(t)  was  derived  in  Chapter  III  and  is  equal  to 


-1 


•5(c9(t))  = dvs  • (s+k2)  • (s  /s”+  as  + k^/i^  b) 


-1 


(52) 
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and 


L (t  1/2)  = r(l/2)//s  = /ir//s. 


(53) 


Therefore, 

i(C2(t))*i(t  l^2)  = d/iKs-H^)  (s/s  + as  + k^/s  + b)  \ (54) 

Using  the  expansion  introduced  in  Chapter  III  (equation  (14)),  equation 
(54)  can  be  expressed  as 


L(C9(t))  • L(t  1/2)  = /7  l 


E. 


(55) 


j=l  /s  • (/i+g  ) 

where  E.  and  g.  are  constants  to  be  determined. 

J J 

Substituting  equation  (55)  into  equation  (51)  and  taking  the  inverse 
Laplace  transform,  one  obtains 


nlow  temp  “ 15/5  • YF(g;]/E)*  (56) 

At  small  t,  the  function  F(g^/t)  can  be  expanded  into  a series  [113] 


F( g./t)  = 1 - — g./tT  + g.  t + 9(t). 

J A 3 


(57) 


Since  £ E.  = 0,  one  obtains 


j=l 


n 


A . 1/2  A . „A, 
‘ + Y1  ' 


low  temp  * Yl/2  * f ' + Yn  • t + 0 (t) 


(58) 


where 


1/2 


p/D  Z (-  — E.g.) 
j=l  /tt  ^ ^ 


(59) 
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and 

A 5 2 

y = p/D  Z (E  gt).  (60) 

j=l  J J 

The  quantity  C2(t)  can  be  expanded  in  a similar  manner  when  t is  small. 
Thus , 


^2  1/2  ^2  ^2 
C2(t)low  temp  ~ Yl/2  + Y1  ' C + 9 (t) 


(61) 


where 


and 


C2  5 2 
Yl/2  " 1-1  ~/7  BjSj} 


C2  5 

\ ‘ * (VjK 

J=1  J J 


(62) 


(63) 


Combining  equations  (58)  and  (61)  with  C^t),  we  obtain  the  alkali  Auger 
intensity  expression  at  low  temperature  and  small  t as 


, .Auger  r_  -kot  . 0 1/2  , 0 , A0.  .. 

X^low  temp  ~ a*-C01  6 + Y1/2C  + Y1  ’ C + 0 


(64) 


where 


0 A ^ C2 

Yl/2  Yl/2  + Yl/2 


(65) 


and 


Yi  - + YX 


(66) 
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Discussion 


A Model  Calculation 

We  have  derived  the  theoretical  expression  for  an  alkali  Auger 
signal  intensity.  Although,  the  reduced  form  is  only  applicable  to  the 
case  of  low  temperature  measurements  for  a small  t.  The  time  dependence 
of  the  intensity  can  be  calculated.  A direct  comparison  of  this  model 
with  the  experimental  data  is  difficult,  however,  a significant  feature 
observed  at  low  temperatures  (the  "incubation  time")  is  theoretically 


predicted. 

The  sodium  Auger  intensity  from  Na20-2SiC>2  has  been  calculated 
using  the  parameters  previously  derived  from  experimental  data  and  liter- 
ature. The  data  are: 

-30  2 

D„  = 10  cm  /sec 

Na 

— 28 

u„  = 10  cm/V  sec 

KNa 

F = 10^  V/cm 

a = 4.0  x 10  ^ cm“  (for  4 keV) 

conv. 

a_cn  = 3.2  x 10  ^ cm^  (for  4 keV) 

ESD 

-30  2 

Note  that  the  diffusion  coefficient  10  cm^/sec  corresponds  to  the 


constant  at  110°K  [104]. 


2 . 


When  a beam  current  density  of  1 mA/cm  is  chosen  for  calculation, 

coefficients  appearing  in  equation  (64)  are  calculated  from  equations 

A 

(65)  and  (66)  and  the  reduced  expressions  of  C^(t)  and  C2(t)  + n (t)  for 
small  t are 


Cl(t)  = CQ1  exp (-0 . 025  t) 


(67) 
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and 

C2(t)  + nA(t)  = CQ  (0.031  t1/2  + 0.010  t + 0(t)).  (68) 

The  alkali  (in  this  case,  sodium)  Auger  signal  intensity  is  proportional 
to  the  sum  of  equations  (67)  and  (68)  as  expressed  in  equation  (37).  It 
is  obvious  that  C^(t)  decreases  exponentially  with  time  while  C2(t)  + 
nA(t)  increases  with  time.  Calculated  curves  from  equations  (67)  and  (68) 

A 

are  plotted  in  Figure  29.  The  quantity  C2(t)  + n (t)  is  only  plotted 

up  to  15  seconds,  since  the  expression  is  valid  only  for  a small  t.  For 

a large  t,  we  have  previously  shown  that  C2(t)  decreases  after  initial 

A 

rise  (Figure  22).  Although  the  expression  for  n (t)  is  complicated,  it 

A 

contains  C2(t)  explicitly  in  the  integrand,  therefore,  n (t)  should  de- 

crease,  similar  to  C9(t)  when  t becomes  large.  Thus,  C2(t)  + n (t)  is 

predicted  to  be  a monotonic  decreasing  function  with  time  for  a large  t. 

A 

However,  from  equation  (68),  C2(t)  + n (t)  increases  with  time  for  a 

small  t.  The  overall  expression  for  the  sodium  Auger  signal  intensity 

is  thus  a sum  of  the  two  curves.  This  is  also  shown  in  Figure  29.  The 

rates  of  decrease  of  equation  (67)  and  increase  of  equation  (68)  are 

approximately  equal  up  to  - 10  seconds.  Therefore,  the  Auger  signal 

intensity  does  not  vary  significantly  for  t < 10  seconds  (equal  to  an 

-2  2 

electron  dose  of  < 1 x 10  coulombs/cm  ).  This  is  recognized  experi- 
mentally and  is  called  an  "incubation  time."  After  this  initial  period 
of  stability,  the  intensity  decays  with  time.  An  incubation  time  is 
calculated  based  on  the  model  involving  the  conversion  process,  the 
electron  stimulated  desorption  and  the  field-enhanced  diffusion.  The 
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Figure  29.  Calculated  sodium  Auger  signal  decay  from  soda-silica  glass. 
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rate  of  conversion  and  desorption  are  controlled  primarily  by  the  beam 
energy  and  the  current  density  but  not  by  the  temperature.  The  rate  of 
diffusion  is  strongly  dependent  upon  the  temperature.  When  the  tempera- 
ture is  high,  the  diffusion  constant  and  the  mobility  are  large  and 
diffusion  is  enhanced.  Consequently,  the  relative  magnitudes  of  C2(t) 

A. 

become  low  and  n (t)  decreases  correspondingly.  This  means  that  the 

A 

rates  of  increase  in  C,,(t)  + n (t)  at  room  temperature  become  smaller 
than  those  calculated  at  liquid  nitrogen  temperature.  Therefore,  the 

rate  of  decrease  in  C^(t)  is  not  equivalent  to  the  rate  of  increase  in 

A A 

C2(t)  + n (t).  The  resultant  curve,  (^(t)  + C2(t)  + n (t),  shows  a mono- 
tonic decrease  and  does  not  have  an  initial  period  at  room  temperature. 
The  alkali  Auger  signal  only  decreases  approximately  exponentially. 

This  is  schematically  illustrated  in  Figure  30. 

When  the  beam  current  density  becomes  high  at  liquid  nitrogen 

temperature,  the  time  response  of  all  the  quantities,  C^t),  C2(t)  and 
A 

n (t)  becomes  shorter  since  the  diffusion  is  enhanced  and  it  limits 

increase  in  ESD  rate.  Therefore,  the  rate  of  initial  rise  in  C2(t)  + 

A 

n (t)  may  not  be  as  large  as  the  rate  of  decrease  in  C^(t).  In  this 
case  the  incubation  time  may  be  shortened  or  may  not  be  observed.  This 
is  consistent  with  the  experimental  observation  shown  in  Figure  5 (from 
[39]).  As  the  beam  current  density  was  increased,  the  incubation  time 
was  shortened.  At  high  current  density  a measurable  incubation  time  was 
not  observed,  but  the  Auger  signal  decayed  approximately  exponentially. 

As  expected,  at  room  temperature  the  incubation  time  was  not  observed. 


Alkali  Auger  Signal  Intensity 
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Figure  30. 


Schematical  illustration  of  relations  for  C^(t),  C2(t) 
and  n^(t)  at  low  temperature  and  room  temperature. 
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Comparison  of  Experimental  Results  with  Model 

The  experimental  sodium  Auger  signal  changes  from  Na20-2Si02  glass 
are  shown  in  Figure  31,  where  curves  (a)  and  (b)  correspond  to  the  measure- 
ments at  110° K and  300° K,  respectively.  A 3 keV  electron  beam  with  1 x 10 
2 

A/ cm  was  used  for  excitation. 

No  incubation  time  was  observed  at  room  temperature.  In  this  case, 
the  decay  curve  is  divided  into  three  parts  (this  is  clearly  seen  when 
the  decay  curve  is  plotted  on  a semi-log  scale) . In  the  first  decay  per- 
iod (first  several  seconds),  a relatively  slow  decay  rate  was  observed. 
Recalling  equations  (37)  and  (64),  this  is  a period  during  which  C2(t)  + 

A 

n (t)  increases,  resulting  in  a very  slow  overall  decay  rate.  However, 

A 

the  quantity  C2(t)  + n (t)  decreases  subsequently  and  the  decay  rate 

becomes  more  rapid  in  the  middle  part.  In  the  third  period  there  is 

almost  no  contribution  from  C2(t)  + n (t) , since  it  asymptotically  appro — 

aches  zero-.  Thus,  the  decay  in  the  third  period  is  controlled  only  by 

changes  in  C^(t).  These  effects  are  illustrated  in  Figure  30. 

Upon  cooling  the  sample  to  110 °K,  an  incubation  time  of  15  seconds 

was  observed  during  which  the  sodium  signal  remained  unchanged.  As  shown 

A 

in  the  model  calculation,  the  increase  in  C2(t)  + n (t)  is  balanced  by 
the  decrease  in  C^(t)  in  this  period.  The  signal  starts  to  decay  subse- 
quently. When  the  beam  current  density  was  changed,  the  length  of  incu- 
bation time  and  decay  rate  varied  correspondingly.  Plotting  the  normalized 
sodium  Auger  signal  versus  electron  dose  (product  of  beam  current  density 

and  time)  in  Figure  32,  all  the  data  essentially  form  a single  line.  The 

-2  2 

sodium  signals  are  unchanged  up  to  1.5  x 10  coulomb/cm  . Beyond  this 


Sodium  Auger  Peak  Height  (a.  u ) 
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Figure  31.  Normalized  sodium  Auger  peak  height  from  soda-disilicate 
glass  as  a function  of  electron  beam  bombardment  time  at 
liquid  nitrogen  and  room  temperatures.  A 3 keV,  1 mA/cm^ 
electron  beam  was  used. 


Normalized  Sodium  Auger  Peak  Height 


96 


Figure  32. 


Normalized  sodium  Auger  peak  height  from  a soda-disili- 
cate glass  as  a function  of  electron  dose  at  liquid 
nitrogen  temperature  (110 °K) . 
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doseage,  the  signal  started  to  decrease.  In  the  model  calculation,  a 

-2  2 

similar  dosage  (~  1 x 10  coulomb/cm  ) is  calculated  for  the  incubation 

period.  Recently,  Pantano  and  Madey  [123]  have  compiled  electron  beam 

damage  data  reported  in  the  literature.  These  are  reproduced  in  Table  2. 

where  the  critical  doses,  D^,  are  those  reported  or  calculated  on  the 

basis  of  some  change  in  the  Auger  spectra  due  to  electron  beam  induced 

damage.  In  many  of  the  references,  the  dose  estimates  are  not  precise, 

so  the  data  are  only  accurate  to  an  order  of  magnitude.  Nevertheless, 

it  is  of  interest  to  note  that  many  organic  and  inorganic  solids  exhibit 

-3  -2  2 

critical  doses  of  10  to  10  coulomb/cm  . This  range  of  values  is  con- 
sistent with  the  dose  which  corresponds  to  the  sodium  Auger  signal  incu- 
bation period  in  Na20-2Si02  glass. 

The  length  of  incubation  time  is  very  dependent  upon  the  glass  compo- 
sition. We  have  observed  an  incubation  time  of  ~ 10  seconds  from  soda- 

2 

disilicate  glass  with  beam  conditions  of  3 keV,  1 mA/cm  . An  incubation 
time  of  more  than  50  seconds  was  reported  by  Pantano  [39]  (Figure  5)  from 
soda-lime-silica  glass  under  nearly  identical  experimental  conditions. 

This  large  difference  in  incubation  time  is  attributed  to  the  structural 
differences  between  the  two  glasses.  In  the  soda-lime-silica  glass,  CaO 
is  substituted  for  Si02  of  a soda-silica  glass.  Sodium  ions  are  then 
bound  more  tightly  in  the  network  [104] and  the  diffusion  constant  of 
sodium  is  much  smaller  [104] . Frischart  [104]  has  compared  the  activation 
energies  for  sodium  self-diffusion  in  soda-silica  and  soda-lime-silica 
glasses.  The  values  were  20  kcal/mole  and  15  kcal/mole  in  the  soda-silica 
glasses  containing  Na20  of  15  and  33  mole  %,  respectively.  Activation 


98 


Table  2.  Critical  Electron  Dose  for  Electron  Bombardment 
Damage  of  Molecular  Solids  and  Films  (from  [123]) 


Material 

Energy 

Dc(coul/ cm^) 

Method 

Elements 

Monitored 

Ref . 

Si3N4 

2 keV 

stable 

AES 

0, Si,N 

[124] 

A1203 

5 keV 

10 

AES 

A1 , 0 

[10] 

Si02 

2 keV 

0.6 

AES 

Si,0 

[124] 

Li2W04 

~1  keV 

0.05 

XPS 

W 

[19] 

NaF,  LiF 

.1  keV 

0.06 

LEED 

[14] 

LiN03,LiS04 

1 keV 

0.05 

XPS, AES 

0,  S ,N 

[19] 

KC1 

1.5  keV 

0.03 

AES 

Cl 

[14] 

Te02 

2 keV 

0.02 

AES 

0 

[124] 

h2o(f) 

1.5  keV 

0.01 

AES 

0 

[20] 

native  oxide 

5 keV 

0.002 

AES 

[9] 

Na3AlF6 

3 keV 

F 

[126] 

CH3OH(F) 

1.5  keV 

2.5  x 10-4 

AES 

C,0 

[20] 
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energies  of  more  than  20  kcal/mole  were  reported  for  soda-lime-silica 
glasses  containing  equal  concentrations  of  ^£0.  This  indicates  that  the 
rates  of  conversion  and  field-enhanced  diffusion  of  sodium  in  the  soda- 
lime-silica  glass  is  much  slower  as  compared  with  the  soda-silica  glass. 

The  sodium  ESD  cross  section  from  soda-lime-silica  glasses  is  expected 
to  be  similar  to  one  from  soda-silica  glasses,  since  the  same  principle 
of  intraatomic  and  interatomic  Auger  de-excitation  can  be  applied.  There- 
fore, the  difference  in  incubation  time  is  due  to  differences  in  the  rates 
of  conversion  and  field-enhanced  diffusion.  Thus,  longer  incubation  times 
are  expected  in  soda-lime-silica  glasses. 

Summary 

1.  A general  mathematical  description  of  the  alkali  signal  decay  has  been 
developed.  It  was  applied  to  low  temperature  measurements  and  the 
origin  of  the  incubation  time  became  apparent. 

2.  Calculations  were  made  of  the  sodium  Auger  signal  measurements  from 
soda-silica  glass  at  110°K.  It  was  demonstrated  that  the  sodium 

_ O 

Auger  signal  was  nearly  constant  for  an  electron  dose  of  - 1 x 10 
2 

coulomb /cm  and  this  was  verified  by  the  experiments. 

3.  Many  organic  and  inorganic  solids  exhibit  electron  damage  (some  change 

— ^ 

in  the  Auger  spectra  due  to  electron  bombardment)  within  doses  of  10 
-2  2 

to  10  coulomb/cm  (from  [123]).  Therefore,  structural  interpretation 
of  the  surface  composition  requires  careful  validation. 

4.  The  length  of  incubation  time  varies  with  the  glass  composition.  Dif- 
ferences in  the  incubation  time  are  attributed  to  structural  differences ; 
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the  higher  the  bonding  strength  and  the  activation  energy  for  alkali 
diffusion,  the  longer  the  incubation  time  is  expected  to  be. 


CHAPTER  V 

APPLICATION  OF  ALKALI  AUGER  SIGNAL  DECAY  PROCESS 
FOR  MATERIAL  CHARACTERIZATION 

In  the  previous  chapters,  the  effects  of  an  electron  beam  on  glass 
during  analysis  by  Auger  electron  spectroscopy  have  been  discussed  and 
have  been  manifested  by  the  alkali  Auger  signal  decay  during  electron 
bombardment.  Non-recognition  of  the  effects  can  seriously  influence  not 
only  quantitative  analysis  but  even  qualitative  analysis. 

However,  these  effects  are  sensitive  to  the  composition  and  structure 
of  the  surface  and  careful  interpretation  of  the  effects  can  provide  fur- 
ther information  about  the  surface  properties.  This  will  be  illustrated 
in  this  chapter,  where  Auger  electron  spectroscopy  will  be  applied  to 
investigation  of  microstructural  characterization  of  glass-ceramics. 

Glass-ceramics  are  recognized  as  extremely  important  engineering 
materials  with  numerous  applications  [127,128].  These  materials  are 
typically  prepared  from  preformed  glass  castings  via  a specific  nucleation- 
crystallization  treatment.  The  controlled  nucleation-crystallization  of 
glass  provides  a microstructure  of  uniform  small  crystals  surrounded  by 
a glassy  matrix.  The  resulting  ceramed  articles  usually  have  superior 
mechanical  and  chemical  properties  as  compared  to  glass  with  the  same 
composition  [127].  The  ceraming  process  produces  local  compositional  and 
structural  alterations  in  the  materials.  These  are  often  difficult  to 
evaluate  with  conventional  x-ray  and  optical  techniques  because  of  the 
lack  of  good  spatial  resolution. 
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In  this  chapter,  it  will  be  shown  that  the  alkali  Auger  signal  decay 
process  may  be  used  to  study  the  microstructure  of  glass-ceramics.  The 
rates  of  alkali  Auger  signal  decay  are  influenced  by  the  presence  of  the 
crystalline  phase.  The  binary  Li20-2Si02  system  was  chosen  in  the  present 
study  because  the  composition  of  the  crystal  is  the  same  as  that  of  the 
bulk  glass  for  this  stoichiometric  material. 

Experimental  Procedure 

Lithia-disilicate  glass  (33  mole  % L^O  and  67  mole  % Si02)  was 
prepared  from  reagent  grade  Id^CO^  and  5 pm  silica  in  the  manner  described 
earlier  in  Chapter  II.  Cast  glasses  were  annealed  at  230° C for  4 hours. 

They  were  nucleated  by  heating  for  24  hours  at  475°C  and  partially  crystal- 
lized under  the  conditions  shown  in  Table  3.  The  respective  volume  fractions 
of  crystallization  of  were  - 30,-  50  and  more  than  70%,  as  inferred 

from  the  data  of  Freiman  and  Hench  [129].  Scanning  electron  microscopy 
was  used  to  observe  fracture  surfaces  etched  with  dilute  HF-.  After  etching, 
the  specimens  were  thoroughly  rinsed,  dried  and  plated  with  gold-palladium 
by  vacuum  deposition. 

The  Auger  measurements  were  made  in  a manner  similar  to  those  described 
in  Chapter  II.  A Physical  Electronics  Thin  Film  Analyzer  ^ was  used.  All 
the  spectra  were  obtained  with  a 1 keV  electron  beam.  The  specimens  were 
cooled  to  - 110°K  by  feeding  liquid  nitrogen  into  a carousel.  The  specimen 


t Cambridge  Stereoscan,  Kent  Cambridge  Scientific,  Inc.,  Morton  Grove,  IL 
it  CMA  10-150,  Perkin  Elmer,  Physical  Electronics  Division,  Eden  Prairie,  MN 
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Table  3.  Specimen  Preparation  Data 


Annealing 

Nucleation 

Crystallization 

Temp(°C) 

Time (h) 

Temp(°C) 

Time(h) 

Temp(°C)  Time (h) 

0% 

glass 

230 

4 

— 

— 

— 

30% 

glass-ceramic 

230 

4 

475 

24 

550  3.3 

50% 

glass-ceramic 

230 

4 

475 

24 

550  8.3 

>70% 



glass-ceramic 

230 

4 

475 

24 

550  24 

104 


temperature  could  be  varied  between  110°K  and  room  temperature  by  simul- 
taneous cooling  with  liquid  nitrogen  and  heating  with  an  electric  heater 
installed  in  the  carousel.  The  temperature  was  measured  with  a Chromel- 
Alumel  thermocouple'  sandwiched  between  two  glass  slides  and  mounted  on 
the  carousel. 


Results 

Scanning  electron  micrographs  of  etched  fracture  surfaces  are  shown 
in  Figure  33.  The  formation  of  rod-like  crystals  reported  by  Freiman  and 
Hench  [129]  is  evident.  The  glass-ceramics  containing  50%  crystal  (b) 
indicates  the  crystal  growth  occurs  to  impingement.  It  appears  that  the 
residual  matrices  apparently  crystallize  in  a dendritic  fashion.  This 
secondary  crystallization  apparently  produces  small  crystallite  within 
the  interstices  of  the  impinged  spherulites  [130].  This  is  no  longer 
evident  in  the  glass  ceramics  containing  more  than  70%  crystal  (a) . The 
glass  phase  apparently  exists  in  continuous  crystal  phases. 

A plot  of  the  lithium  Auger  peak  height  variation  as  a function  of 
beam  bombardment  time  for  the  uncrystallized  glass  is  shown  in  Figure  34. 

The  data  were  taken  at  - 110°K.  As  expected,  the  decay  rates  were  strongly 
dependent  on  the  beam  current  densities  and  faster  decay  rates  were  observed 
with  higher  beam  current  densities.  Incubation  times  were  observed  for 
beam  current  densities  of  5 1 x 10  A/cm  . For  example,  at  7 x 10  A/cm  , 

the  incubation  time  was  about  20  seconds  (equal  to  an  electron  does  of 

-2  2 
~ 1.4  x 10  coulomb/cm  ). 


t Hoskins  Manufacturing  Company,  Detroit,  MI 
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(a)  more  than  70% 


(b)~50% 


(mag.  xlOOO) 


Figure  33.  Scanning  electron  micrographs  of  etched  fracture  surface 
of  lithi-disilicate  glass  containing:  (a)  more  than  70% 

crystallization  and  (b)  approximately  50%  crystallization. 


Normalized  Lj  Auger  Peak  Height 
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Beam  Exposure  Time(min.) 


Figure  34.  Lithium  Auger  peak  height  as  a function  of  beam  exposure 
time  and  beam  current  density  for  lithia-disilicate  glass 
at  liquid  nitrogen  temperature  (110°K). 
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The  lithium  Auger  signal  versus  irradition  time  from  glass  ceramics 
with  30%  (by  volume)  crystallization  is  shown  in  Figure  35.  Incubation 
times  similar  to  the  uncrystallized  glass  were  observed.  As  the 
beam  current  density  was  increased,  the  decay  rate  was  accelerated  as 
expected,  but  there  were  additional  features  in  the  decay  which  were  not 
observed  in  the  uncrystallized  glass.  Three  regions  of  behavior  in  the 
lithium  Auger  signal  decay  were  observed,  as  shown  in  Figure  35.  In 
region  I,  the  signal  for  the  glass-ceramic  decayed  at  approximately  the 
same  rate  as  the  glass  for  identical  beam  conditions.  The  rates  of  signal 
decays,  however,  became  considerably  slower  in  glass-ceramics  for  all  beam 
current  densities,  as  indicated  in  region  II.  For  a current  density  of 

— O O 

1.25  x 10  A/ cm  the  signal  was  nearly  independent  of  beam  irradiation 
time.  In  region  III,  the  signal  again  decayed,  but  the  rate  was  less 
than  for  glass.  The  extent  of  each  region  was  strongly  dependent  upon 
the  beam  current  density. 

Similar  results  were  obtained  for  the  glass-ceramics  with  higher 
volume  fraction  crystallinity.  The  decay  curves  for  four  volume  fractions 
of  crystallization  (0%,  -30%,  -50%  and  more  than  79%  crystalline)  are 
shown  in  Figure  36.  These  data  were  obtained  with  a current  density  of 
1.25  x 10  A/cm  at  - 110°K.  For  each  glass-ceramic  specimen,  the 
lithium  Auger  signal  decayed  rapidly  in  region  I,  followed  by  a nearly 
constant  signal  which  was  independent  of  the  beam  irradiation  time  (region 
II).  Decay  was  again  observed  in  region  III. 

The  temperature  dependence  in  the  lithium  Auger  signal  decay  from 
the  '30%  crystallized  specimen  is  shown  in  Figure  37.  A 1 keV,  1.25  x 10 


Normalized  L|  Auger  Peak  Height 
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Lithium  Auger  peak  height  as  a function  of  beam  exposure 
time  and  beam  current  density  for  lithia-disilicate  glass 
containing  30%  crystallization  at  liquid  nitrogen  temper- 
ature. 


Figure  35. 
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Beam  Exposure  Time  (min.) 


Figure  36.  Lithium  Auger  peak  height  as  a function  of  beam  exposure 
time  for  various  volume  fractions  of  crystallization. 
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Figure  37.  Lithium  Auger  peak  height  as  a function  of  substrate 

temperature  for  33L  glass  containing  30%  crystallization. 


Ill 


2 

A/ cm  electron  beam  was  used.  A length  of  region  II  progressively 
decreased  as  the  specimen  temperature  increased. 

Discussion 

The  lithium  Auger  signal  decays  from  glass  and  from  glass-ceramics 

were  considerably  different.  It  was  shown  that  the  decay  rates  from 

glass-ceramics  were  nearly  identical  to  those  from  glass  during  early 

electron  bombardment  but  became  significantly  different  upon  further 

bombardment.  These  were  clearly  shown  by  comparing  the  Auger  decays 

from  uncrystallized  glass  (Figure  34)  with  those  from  30%  crystallized 

glass  (Figure  35).  In  both  specimens,  nearly  identical  decay  rates  were 

observed  in  the  first  60  seconds  (region  I) . These  data  suggest  that 

the  lithium  ions  in  the  glass  phase  of  the  glass-ceramics  are  responsible 

for  the  Auger  signal  loss.  After  the  lithium  ions  in  the  glass  phase 

were  depleted,  the  Auger  signals  in  region  II  are  due  primarily  to  the 

lithium  ions  in  the  crystal  phase.  These  features  are  further  substan- 

_3  2 

tiated  in  Figure  36,  when  a beam  current  density  of  1.25  x 10  A/cm 

was  chosen.  The  Auger  signal  remained  relatively  constant  in  region  II 

and  its  normalized  value  corresponds  to  the  percent  of  crystallization. 

We  will  define  this  period  as  a "ceramic  incubation  time,"  similar  to 

an  incubation  time  observed  in  glass. 

The  existence  of  a relatively  long  "ceramic  incubation  time"  for 

the  crystalline  phase  is  reasonable  based  upon  its  properties.  For 

example,  (1)  the  diffusion  coefficient  of  Li  in  the  crystalline  phase 
2 3 

is  at  least  10  to  10  times  lower  than  that  in  the  glass  of  the  same 
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composition  [131],  (2)  the  crystalline  phase  has  a higher  thermal  conduc- 
tivity than  the  glass  phase  of  the  same  composition  [131,132]  and  (3)  the 
bond  between  the  Li  ion  and  the  non-bridging  oxygen  is  stronger  in  the 
crystalline  phase  than  in  the  glass  phase  [131].  From  these  facts,  the 
ceramic  incubation  time  is  interpreted  from  the  discussion  described  in 
Chapter  IV  in  a qualitative  sense.  For  the  crystalline  phase,  because 
of  the  higher  bond  strength  between  Li  and  0,  the  rate  of  Si^+0^  Li‘'"+ 
defect  formation,  or  the  conversion  rate,  k2»  is  low.  Thus,  the  exponen- 
tial decay  of  C^(t)  in  equation  (37)  is  slow.  Local  heating  due  to  the 
electron  beam  is  expected  to  be  smaller  for  the  same  power  input  because 
of  the  higher  thermal  conductivity.  The  field-enhanced  diffusion  is  also 

smaller  because  of  the  lower  diffusion  coefficient.  Therefore,  the 

A 

increase  in  C9(t)  + n (t)  is  extended  to  a longer  bombardment  time  and 
it  is  balanced  with  the  decrease  in  C^(t).  Consequently,  this  results 
in  a longer  incubation  time  in  the  crystalline  phase.  This  period 
becomes  progressively  shorter  when  the  beam  current  density  is  increased 
(Figure  35).  In  this  case,  the  area  being  analyzed  in  the  specimen  is 
subjected  to  higher  temperature  due  to  the  beam  heating  and  the  mobility 
of  Li  increases.  In  addition,  the  electron  stimulated  desorption  becomes 
important  as  the  beam  current  density  is  increased  (in  the  case  of  Na, 
see  Figure  26).  Therefore,  the  decrease  in  ceramic  incubation  time  is 
attributed  to  an  increase  in  the  rate  of  diffusion  and  ESD  of  Li. 

Although  the  existence  of  relatively  long  incubation  times  was 
reasonably  explained,  the  ceramic  incubation  time  was  shortened  as  the 
volume  fraction  of  crystallization  increased  (Figure  36) . This  is  proposed 
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to  result  from  the  power  dissipation  in  the  partially  crystallized  speci- 
mens. Although  the  thermal  conductivity  in  the  crystalline  phase  is 
higher,  an  increase  in  the  crystallized  volume  fraction  introduces  poro- 
sity (due  to  thermal  stress  at  the  phase  boundary) , as  shown  in  the 
scanning  electron  micrographs  (Figure  33).  Thereby,  the  overall  thermal 
conductivity  of  the  composite  decreases  due  to  increased  phonon  scattering 
and  a lower  conduction  heat  flux.  Therefore,  the  local  temperature  from 
the  electron  beam  would  increase  as  the  volume  fraction  of  crystallization 
increased,  with  other  factors  remaining  constant.  Besides  a higher  local 
temperature,  the  deposited  charges  (electrons)  tend  to  accumulate  at  phase 
boundaries  of  individual  crystallites  [115].  The  amount  of  charge  accumu- 
lation is  thereby  increased  and  a higher  electric  field  may  result. 

Higher  local  temperatures  and  higher  electric  fields  increase  Li  mobility 
and  the  diffusion  driving  force,  which  results  in  an  increased  field-enhanced 
diffusion  rate.  The  ceramic  incubation  time  observed  in  specimens  with  a 
high  volume  fraction  of  crystallization  is  therefore  decreased.  This 
hypothesis  was  tested  by  looking  at  the  ceramic  incubation  time  as  a 
function  of  the  substrate  temperature,  as  shown  in  Figure  37.  As  the 
temperature  of  the  specimen  increased,  the  local  temperature  correspon- 
dingly increased  for  a constant  beam  power.  Thus,  the  lithium  mobility 
increased  and  the  field-enhanced  diffusion  rate  increased  correspondingly, 
resulting  in  shorter  ceramic  incubation  times  at  higher  temperatures. 

The  discussion  presented  above  shows  that  thermal  conduction  plays 
an  important  role  in  the  rate  of  the  alkali  Auger  signal  decay  in  the 
composite  material.  The  thermal  conductivity  is  strongly  dependent  upon 
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the  specimen  microstructure;  therefore,  careful  interpretation  is  neces- 
sary for  Auger  analysis  of  the  composite  material.  Nevertheless,  the 
present  study  indicates  that  the  electron  beam  effects  can  be  used  to 
provide  a basis  for  Auger  electron  spectroscopic  studies  of  microstructural, 
as  well  as  compositional  properties. 

Summary 

1.  Lithium  Auger  signal  variation  has  been  studied  from  the  lithia- 
disilicate  glasses  containing  L^S^O,.  crystallite  with  a volume 
fraction  of  0%,  -30%,  ~50%  and  more  than  70%. 

2.  The  rates  of  decay  of  the  lithium  Auger  signals  from  crystallized 
glasses  were  influenced  by  the  presence  of  the  crystalline  phase 
and  were  strongly  dependent  on  beam  current  density. 

3.  Upon  choosing  an  appropriate  beam  current  density,  the  decay  curve 

from  crystallized  glasses  was  divided  into  three  regions  with  respect 
to  electron  bombardment  time.  (a)  Region  I:  the  decay  characteristics 

were  identical  to  those  of  uncrystallized  glass  and  were  attributed  to 
the  lithium  Auger  signal  decay  in  the  glass  phase.  (b)  Region  II: 

the  lithium  Auger  signal  remained  constant  and  its  normalized  inten- 
sity corresponded  to  the  percent  of  crystallization.  This  period  was 
named  "ceramic  incubation  time."  (c)  Region  III:  this  lithium  Auger 

signal  decay  was  attributed  to  the  crystalline  phase  of  the  material. 

4.  The  ceramic  incubation  time  was  interpreted  from  the  general  alkali 
Auger  signal  decay  theory  developed  in  Chapter  IV. 
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5.  The  ceramic  incubation  time  decreased  as  the  volume  fraction  of 
crystallization  increased.  This  was  explained  by  the  power  dissi- 
pation and  the  charge  accumulation  in  the  partially  crystallized 
specimens . 

6.  The  present  study  indicates  that  the  electron  beam  effects  can  be 
used  to  provide  a basis  for  Auger  electron  spectroscopic  studies 
of  microstructure  as  well  as  composition. 


CHAPTER  VI 

SUMMARY  AND  CONCLUSIONS 


The  goal  of  the  present  research  was  to  investigate  possible  mechanisms 
of  alkali  Auger  signal  loss  on  glass  surfaces  during  electron  bombardment. 

It  was  concluded  that  the  signal  loss  was  caused  by  the  disappearance  of 
alkali  ions  at  the  surface  due  to  electron  stimulated  desorption  and  elec- 
tric field  enhanced  diffusion.  This  was  modeled  by  considering  the  following 
processes:  An  energetic  electron  will  break  a bond  to  create  an  ion  pair 

of  oxygen  and  alkali  ions  at  the  surface.  The  oxygen  ion  will  then  trap 
a hole  and  the  alkali  ion  will  no  longer  be  bound.  Further  bombardment  by 
electrons  will  cause  an  electronic  excitation  of  oxygen  (de-excitation) 
followed  by  an  escape  of  the  alkali  ion  via  coulombic  repulsion  (electron 
stimulated  desorption).  In  addition,  an  electric  field  is  produced  normal 
to  the  surface  as  a result  of  the  surface  positive  potential  and  the  diffused 
negative  layer  at  a depth  equal  to  the  electron  range.  Under  the  influence 
of  this  electric  field  the  alkali  ions  diffuse  away  from  the  surface  (elec- 
tric field  enhanced  diffusion) . This  is  further  accelerated  when  electron 
beam  heating  causes  an  increase  in  temperature,  resulting  in  higher  alkali 
mobility. 

This  general  model  was  then  formulated  mathematically  and  a closed 
form  of  analytical  solution  was  found.  The  solution  includes  all  the 
experimental  parameters,  such  as  electron  beam  energy,  current  density, 
bombardment  time  and  temperature;  and  the  material  parameters,  such  as 


116 


117 


diffusion  constant  and  mobility.  Relative  contributions  of  desorption 

and  diffusion  to  the  processes  were  calculated  and  it  was  found  that  ESD 

-3  2 

becomes  important  with  beam  current  densities  greater  than  10  A/cm 

but  the  diffusion  limited  further  increase  in  ESD  with  densities  greater 
-2  2 

than  5 x 10  A/cm  . In  a conventional  AES  analysis  of  glass  surfaces, 

both  desorption  and  diffusion  are  occurring  and  are  important. 

An  incubation  time  during  which  the  alkali  Auger  signal  is  constant 

was  demonstrated  by  the  general  model,  when  the  solution  was  applied  to 

the  case  of  low  temperature.  For  soda-silica  glass,  the  sodium  Auger 

-2  2 

signal  was  nearly  consistent  within  a dose  of  ~ 1 x 10  coulomb/cm  and 

this  was  verified  experimentally.  In  this  case,  the  incubation  time  was 

-3  2 

only  several  seconds  when  an  electron  beam  with  ~ 10  A/ cm  range  was 

used.  It  is  emphasized  that  the  analysis  has  to  be  made  within  this  time 
after  the  electron  gun  was  turned  on  and  the  data  taken  after  this  time 
is  not  meaningful.  The  reluctance  to  recognize  the  effect  can  seriously 
influence  the  analysis.  The  length  of  the  incubation  time  is  also  depen- 
dent upon  the  glass  composition.  It  follows  from  theory  that  glasses 
containing  alkali  ions  with  lower  diffusion  constants  and  higher  bonding 
strengths  result  in  longer  incubation  times  (for  example,  soda-lime-silica 
glass).  Thus,  the  theory  developed  in  the  present  work  has  a practical 
importance  for  glass  surface  analysis  by  AES. 

Electric  field  enhanced  diffusion  of  sodium  was  demonstrated  in  soda- 
silica  glass  thin  films  deposited  on  stainless  steel  substrates.  With 
these  samples,  the  amount  of  charge  trapping  was  controlled  and  the  region 
of  the  glass  over  which  the  diffusion  would  occur  was  defined.  Thus,  a 
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mathematical  model  was  developed  and  compared  with  the  experimental  data. 

The  electron  beam  raised  the  temperature  of  the  analysis  site  by  approxi- 
mately 200°C  and  the  field  strength  was  determined  to  be  approximately 
10"*  V/cm.  The  model  was  also  used  to  calculate  the  distribution  of  sodium 
with  depth  in  the  film,  but  the  calculated  distribution  disagreed  with 
experimental  sputter  profile  data.  This  results  from  redistribution  of 
sodium  during  Ar+  bombardment  for  sputter  profiling. 

Sodium  desorption  kinetics  were  studied  from  soda-silica  glass  surfaces 

fractured  in  situ.  It  was  well  explained  by  the  general  model  described 

-19  2 

above.  The  ESD  cross  section  was  evaluated  to  be  approximately  10  cm 
and  decreased  as  the  primary  electron  energy  increased.  A conventional 
ESD  mechanism  (Menzel-Gomer-Redhead  model)  was  not  adequate  to  explain 
the  sodium  desorption  from  glass.  In  this  case,  sodium  desorbs  via  intra- 
atomic  Auger  de-excitation  for  the  primary  electron  energies  of  less  than 
100  eV,  while  interatomic  Auger  de-excitation  is  the  most  likely  mechanism 
above  100  eV.  For  higher  energies  (e.g.,  > 530  eV),  it  is  speculated  that 
the  sodium  desorbs  via  both  intra  and  interatomic  Auger  de-excitation  pro- 
cesses. The  ESD  cross  section  was  correlated  with  these  mechanisms  and 
the  experimentally  determined  ESD  cross  section  was  shown  to  be  very 
reasonable.  Because  of  unique  characteristics  of  the  ESD  mechanism,  it 
is  suggested  that  ESD  can  be  a valence-sensitive  probe  of  surfaces  and 
careful  measurements  of  ESD  threshold  energy  and  desorbed  ion  energy 
distribution  may  provide  further  information  about  the  surface.  Further 
investigation  of  ESD  with  this  point  in  mind  may  prove  to  be  very  useful. 
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Finally,  AES  was  applied  to  the  investigation  of  microstructural 
characterization  of  lithia-silica  glass-ceramics.  The  rates  of  lithium 
Auger  signal  decay  were  influenced  by  the  presence  of  the  crystalline 
phase.  With  proper  experimental  conditions,  the  decays  from  glass  and 
crystalline  phases  were  differentiated  and  the  normalized  lithium  Auger 
peak  height  from  the  crystalline  phase  corresponded  to  the  percent  of 
crystallization.  This  feature  indicates  a basis  for  AES  studies  of  micro- 
structural  as  well  as  compositional  properties  of  materials. 


APPENDIX  A 

SOLUTION  OF  EQUATION  (3) 

Equation  (3)  will  be  solved  for  the  boundary  conditions  defined  by 
equations  (6)  and  (7). 

If  a solution  of  equation  (3)  is  assumed  to  be  equal  to  the  product 
of  a function  Z of  z,  only  distance,  and  by  a function  T of  only  time,  t. 


one  obtains 

T'/T  = (DZ"  - EZ')/Z  = constant  = -t 

(A-l) 

where  E = F • y and  r is  as  yet  unrestricted.  Solutions  for  T and  Z are 
av  J 


T = exp  (-rt) 

(A-2) 

Z = exp (Ez/2D) • (K^  sin  coz  + cos  coz) 

where 

(A-3) 

co  = (4rD  - E2)1//2/2D. 

(A-4) 

Solving  equation  (A-4)  for  r in  terms  of  w,  substituting  this  into 
and  taking  the  product  of  (A-2)  and  (A-3),  one  obtains 

(A-2), 

2 2 2 

c(z,t)  = exp{-(4D  w“+E  )t/4D}*exp(Ez/2D)*  (K^  sin  wz+I^cos  uz)  . 

(A- 5) 

The  substitution  of  (A-3)  in  the  boundary  conditions  leads  to  the  equations 


K a»D  - 1/2  K2E 

(A-6) 

K (m-D  cos  ul  - 1/2E  sin  coil  )-K2  (1/2E  cos  uJl+coD  sin  coil)  = 0. 

(A- 7) 
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In  order  that  these  homogeneous  equations  possess  solutions  and  Kj 
other  than  = 0,  it  is  necessary  that  the  determinant  of  the 

coefficient  vanishes.  That  is 

sin  • (u2D2  + E2/4)  = 0 

so  that 

to  = nnr/Jl  m = 0,1,2...  (A-8) 

or 

co  - ± £E/2D.  (A- 9) 

The  constant  u>  has  one  of  the  values  given  by  (A-8)  or  (A-9).  The  ratio 
of  to  K2  is  given  by  either  equations  (A-7)  or  (A-6) . If  one  chooses, 

K = K E 

1 m 

then 

K0  = K ^nnrD/i. 

2 m 

When  m has  the  value  given  by  (A-9),  r is  zero,  so  that  T is  unity  and 
the  corresponding  special  solution  of  equation  (3)  is 

c = Cq  exp (Ez/D) . 

Thus,  from  (A-5) , the  solution 

00 

c = exp  {-(4D2m2iT2+E22,2)t/4D£2}  • exp(Ez/2D)*  I K^fE  sin  (nnrz/i)  + 
(2mirD/5,)  • cos  (mirz/i  ) ] + Kq  exp  (Ez/D) 


(A-10) 
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satisfies  equation  (3)  and  conditions  (6)  and  (7),  where  and  have 
not  yet  been  determined.  For  t=0,  equation  (A-10)  gives 


c exp(-Ez/2D)  = K exp(Ez/2D)+  £ K [E  sin(m7rz/£)+(2miTD/£) 


• cos(nnrz/£)] 


(A-ll) 


and  it  is  necessary  to  determine  the  constants  K and  K so  that  this 

0 m 

equation  will  be  satisfied.  Let 

u = E sin  to  z + 2Dto  • cos  to  z (A-12) 

m m m m 

and 


uQ  = exp (Ez/2D) 

where 


2 2 2 

u>  = mr /£  m ^ 0 and  to.  = -E  /4A  . 

m 0 

The  function  u^  satifies  the  differential  equation  of 

u'.'  + uk  u . = 0 i = 0 , 1 , 2 . . . 
ill  ’ ’ 

and  the  boundary  conditions 

u ! = E/ 2D  • u . z = 0 and  £ . 

l l 

Hence,  integrating  by  parts,  we  obtain, 


£ 

I 

0 


(u'.'  u.-u'.'  u.)dz  = [u!  u . -u ! u . - /(u!u!-u!u!)dz  = 0. 

1 j 3 i 1JliOnvijjiy 


(A-13) 


(A-14) 
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Prom  equation  (A-14),  it  follows  that 


/ (u'.'u . -u'.'u  ,)dz  = (w2-u)2)  / u.u.dz. 

0 1 J J 1 J 1 o 1 3 


Hence,  the  integral  from  0 to  i of  the  product  u.u. 

i 3 


should  vanish  when 


i t4  j and  only  when  i = j , 


f u2dz  = 1/2 (E2  + 4D2w2),  i^O.  (A-15) 

0 1 1 

Using  these  features,  multiplying  both  sides  of  equation  (A-ll)  by  u^, 
and  integrating  from  0 to  l yields 

a i 

K_  f exp(Ez/D)dz  = / c dz  = c i 

0 0 0 ° 

then 


Kq  = EcQJi/D{exp(EZ,/D)-l}. 

Similarly,  multiplying  both  sides  of  equation  (A-ll)  by  u when  m i 0 

m 

and  integrating  from  0 to  £ , one  obtains 

Km  = 16D2Emr£3.  {4D2m2ir2+E2£2 }"2  • {1  ± exp(-E£/2D)} 


with  the  minus  sign  for  m even,  and  the  plus  sign  for  m odd.  Substituting 
Kq  and  into  equation  (A-10) , the  final  solution  of  equation  (3)  with 
conditions  (6)  and  (7)  is  obtained. 


APPENDIX  B 

SOLUTION  OF  EQUATION  (12) 


In  this  appendix,  equation  (12)  will  be  solved.  It  is  possible  to 
reduce  equation  (12)  to  an  ordinary  diffusion  equation  by  means  of  the 
following  transformations  [133], 


:b(z,t)  = C*(z,t)  • expt^z  - ^-t  ]. 


(B-l) 


The  reduced  form  is  then 

Jb 


3Cfa(z 


3t 


,t)  a2c*(z,t) 

= D 


(B— 2 ) 


3 z 


and  the  boundary  conditions  are 


C*(0,t)  = p cz(t)  exp(^-  t) 


and 


C*(0,0)  = 0 (B-3) 

Introducing  the  Laplace  transformation 
C* ( z , s ) = L (C* (z , t) } 

and 

C2*(z,s)  = I{C2(z,t)}. 


Equation  (B-2)  becomes 
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Hence,  using  the  boundary  conditions  (B-3) , one  obtains 

2 


C*(2,s)  - PC2  -(s  - SHI 


) • exp(-z/s/D)  . 


(B-4) 


The  inverse  Laplace  tranf ormation  of  the  exponential  is 
2 

L ^ exp(-z/s/D)  = <*exp(-  ^-)*  (irt3) 

where 


K = 


2/d- 

Using  the  face  that 


r~  1 f/i  / _ (f  y ) \ 1 

L {C2(S  " ~ 4D  * 


C2(t)  • exp  (■—-  t) 


and  the  convolution  theorem  (Appendix  C)  for  Laplace  transformation,  one 

* 

obtains 

t 2 

C*(z,t)  = p / C2(t')-  exp ( ~ t ) • k • exp(-<2(t-t')  1)  (iT(t-t')3) 

° (B-5) 


APPENDIX  C 
CONVOLUTION  THEOREM 

Theorem  [134] 

If  e S0Cf^(t),  e S0Cf2(t)  are  bounded  in  t > 
t 

L{J  f1(t-t')  • f2(t')dt’}  = 0L-02 

where 

0-^s)  = Lif^t)}  and  02(s)  = i (f  2 (t ) } . 


0 , then 
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APPENDIX  D 

ELECTRON  BEAM  HEATING 


In  order  to  estimate  the  temperature  in  the  film,  the  heat  transport 
shown  in  the  below  is  considered. 


For  an  arbitrary  radius  r_.  the  condition  of  heat  balance  is 
r 

/ [p  (r ' )-w(r ' ) ] 2tt  r'dr'  - s(r)27rrh  = 0 (D-l) 

0 

From  this,  the  temperature  distribution  can  be  calculated  using  w(r)  = 

- a (3T(r)/3Z)  and  S = - a (3T(r)/3r),  where  a and  a are  the  thermal 
g m g m 

conductivities  of  glass  and  metal,  respectively.  Equation  (D-l)  can  be 
solved  approximately  for  small  uniform  electron  beam  (radius  to  < L) . The 
temperature  rise  AT  obtained  in  the  center  of  the  beam  is  given  by 

AT  = pr  /a  / (1  + 0.67  (a  /a  )(h/r  ) (D-2) 

° g g m o 

Assuming  rn  = 500  q,  a = 2.5  W cm  ^ K \ a = 0.01  W • cm  ^ K ^ and  the 
o m g 

-2 

input  power  of  20W  cm  , the  temperature  is  expected  to  increase  approxi- 
mately  150  C at  the  film  thickness  of  3000  A.  At  this  temperature  rise, 
the  sublimation  of  sodium  from  glass  is  negligible. 
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